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ABSTRACT 

We present the ATLAS^^ Hi survey of a volume-limited, complete sample of 166 nearby 
early-type galaxies (ETGs) brighter than Mk= -21.5. The survey is mostly based on data 
taken with the Westerbork Synthesis Radio Telescope, which enables us to detect H I down to 
5 X 10^ - 5 X 10^ M© within the survey volume. 

We detect H I in ~ 40 percent of all ETGs outside the Virgo galaxy cluster and in ~ 10 
percent of all ETGs inside it. This demonstrates that it is common for non-cluster ETGs to host 
H I. The morphology of the detected gas varies in a continuous way from regular, settled H I 
discs and rings to unsettled gas distributions (including tidal- or accretion tails) and systems 
of gas clouds scattered around the galaxy. The majority of the detections consist of H I discs 
or rings (1/4 of all ETGs outside Virgo) so that if H I is detected in an ETG it is most likely 
distributed on a settled configuration. These systems come in two main types: small discs 
(M(H l)< 10^ Mq), which are confined within the stellar body and share the same kinematics 
of the stars; and large discs/rings (M(H l) up to 5 x 10^ M©), which extend to tens of kpc from 
the host galaxy and are in half of the cases kinematically decoupled from the stars. 

Neutral hydrogen seems to provide material for star formation in ETGs. Galaxies con- 
taining H I within ~ \Rq exhibit signatures of on-going star formation in ~ 70 percent of the 
cases, ~ 5 times more frequently than galaxies without central H I. The ISM in the centre of 
these galaxies is dominated by molecular gas. In ETGs with a small gas disc the conversion 
of H I into H2 is as efficient as in spirals. 

The ETG H I mass function is characterised by M* ~ 2x 10^ M© and by a slope a 0.7. 

Compared to spirals, ETGs host much less H I as a family. However, a significant fraction of 
all ETGs are as H l-rich as spiral galaxies. The main difference between ETGs and spirals is 
that the former lack the high-column-density H I typical of the bright stellar disc of the latter. 

The ETG H I properties vary with environment density in a more continuous way than 
suggested by the known Virgo vs. non- Virgo dichotomy. We find an envelope of decreasing 
M(H l) and M(H i)/Lk with increasing environment density. The gas-richest galaxies live in 
the poorest environments (as found also with CO observations), where the detection rate of 
star-formation signatures is higher. Galaxies in the centre of Virgo have the lowest H I content, 
while galaxies at the outskirts of Virgo represent a transition region and can contain significant 
amounts of H I, indicating that at least a fraction of them has joined the cluster only recently 
after pre-processing in groups. Finally, we find an H I morphology-density relation such that at 
low environment density (measured on a local scale) the detected H I is mostly distributed on 
large, regular discs and rings, while more disturbed H I morphologies dominate environment 
densities typical of rich groups. This confirms the importance of processes occurring on a 
galaxy-group scale for the evolution of ETGs. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: ISM - galaxies: evolution 



1 INTRODUCTION 

In the classification scheme proposed by |Hubble| ( |1936| ) galaxies 
are arranged on a sequence going from ellipticals to lenticulars and, 
from these, to spirals of progressively later type (Sa to Sc). Ellip- 
ticals are dominated by a stellar bulge while the spiral sequence 
is essentially one of decreasing bulge-to-disc ratio. The interme- 
diate position of lenticular galaxies in this scheme has lead to the 
common idea that all early-type galaxies (ellipticals and lenticulars; 
hereafter ETGs) have higher bulge-to-disc ratio than spirals. 

The other difference between ETGs and spirals is that the for- 
mer lack the blue spiral arms typical of the latter (Hubble|1926| ). 
It was early recognised that this corresponds to a lack of star for- 
mation in ETGs, leading to the simplified picture that their stellar 
populations are uniformly old. 

In contrast with this traditional view, ETGs exhibit a large va- 
riety of shapes and some authors suggest that their bulge-to-disc 
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ratio can in fact be as low as that of Sc, disc-dominated spirals 
fSpitzer & Baade|p^5T] [Sandage, Freeman & Stokes||1970[ [van 
den Bergh 1 976| ). Furthermore, following early insights by |Gunn, 
Stryker & Tinsley ( 1981 ) and Rabin] ( |1982| |, it is now estabhshed 
that a large fraction of ETGs are forming small amounts of stars or 



have done so in their recent past (e.g., |Gonzalez| 1993 


JTrager et al.| 


|2000||Yi et al.|2005||Kaviraj et al.|2007||Thomas et a] 


.|2010|). 



Support to these ideas comes from recent studies of nearby 
ETGs using integral-field spectroscopy { |de Zeeuw et al.|[2002| ). 
These show that most ETGs host a rotating, kinematically cold 
component ( [Krajnovic et al.|2008 ) whose stars are usually younger 
and more metal-rich than those in the bulge ( [Kuntschner et al.j 
I20T0I ). 

These results are placed on a firm statistical ground by 
ATLAS^^, a multi- wavelength study of a volume-limited sample of 
260 morphologically selected ETGs ( [Cappellari et al.|201 1al. here- 
after Paper I). We find that as many as 80 percent of all ETGs in the 
nearby Universe consist of nearly axisymmetric, fast rotating stellar 
systems (Krajno vic et al.|2()TT] herafter Paper II; [Emsellem et al.j 
|2011| hereafter Paper III), most of which resemble spiral galaxies 
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with the arms removed ( [Cappellari et al.|[2011b| hereafter Paper 
VII). 

The presence of a substantial disc and the occurrence of star 
formation in ETGs imply that cold gas has played an important role 
in their evolution. Indeed, Khochfar et al. (2011 ) suggest that most 
ETGs grow a stellar disc following gas cooling. In this respect, two 
fundamental lines of research are the direct observation of neutral 
hydrogen gas (Hi) and molecular gas (H2, observed through the 
radiation emitted by CO molecules). In spirals, H I is an essential 
constituent of the disc, and is the material from which H2 and sub- 
sequently new stars form. Understanding the H I and H2 proper- 
ties of ETGs is therefore crucial to investigate the origin of their 
structure and star formation history, how they continue to evolve at 
z = 0, and why they are so different from spiral galaxies. 



etaLLl986, Lake, Schommer & van Gorkom 1987, Shostak 1987, 
Kim et al. 1988 , van Driel & van Woerden 1991, Schiminovich] 



Results from the ATLAS^^ CO survey are presented in ( Young 
|et aL|20l7] hereafter Paper IV). Here we present an H I survey of 
166 nearby ETGs belonging to the ATLAS^^ sample. Thanks to its 
unprecedented combination of sample size, sample selection, depth 
and resolution of the H I data, and availability of multi- wavelength 
data this survey represents a significant improvement over previous 
studies. 

The study of H I in ETGs dates back to the end of the 1960's. 
Early work already showed that ETGs have lower M(H I) /Lb than 
spirals, consistent with their redder colour (e.g ., [G ouguenheim 



1969[|Bottinel li et al. 119701 IB alkowski et al.|1972[ [Gariaghe r^l972( 
Bottinelli, G ouguenheim & Heidmann| |1973| [ Gallagher, F aber &I 



Balick 1975, Balick, Faber & Gallagher 1976, Bieging & Bier- 
mann 1977, Knapp et al. 1977 , Knapp, Kerr & WiUiams 1978, 
Kna pp, Kerr & Hender s on| 1 979 [ [B iermann, Clarke & Frick e|1979[ 
^mm&Sa lpeter 1979;"Sanders 1980, Reif et al. 198 2j. Based on 
these data, [Knapp, Turner & Cunniffe ( 1985 ) and Ward le & Knapp | 
(1986) analysed the Hi content of ~ 150 ellipticals and ~ 300 
lenticulars, respectively. They detected Hi in about 15 percent of 
all ellipticals and 25 percent of all lenticulars for a typical M(H I) 
detection limit of a few times 10^ M©. They also found a lack of 
correlation between M(H I) and Lb and interpreted it in terms of 
gas external origin. 

Early studies were mostly carried out with pointed single-dish 
observations aimed at determining global quantities like integrated 
H I mass and velocity width of ETGs. Recent progress in this kind 
of analysis has been made possible by single-dish blind surveys 
of large areas of the sky such as HIPASS ( [Barnes et al.|200l| see 
[Sadler, Oosterloo & Morganti|2002|for an H I study of ETGs based 
on these data) and ALFALFA ( [Giovanelli et al.|2005 |. In particu- 
lar, the latter pushes the M(Hl) detection limit below 10^ M© for 
galaxies within a few tens of Mpc from us, allowing a significant 
increase in the number of detected ETGs. ALFALFA data show 
that the ETG H I detection rate is a strong function of environment 
density, from just a few percent inside the Virgo galaxy clus ter ([di] 
Serego Alighieri et al.'2007^ to about 40 percent outside it ( Grossi 
et al. 2009, see also jChamaraux, Balkowski & Fontanelli 1986). 



This result fits with the idea that H I is easily stripped from galax- 
ies inside clusters (e.g., [Giovanelli & Haynes|1983| ), and with the 
fact that recent episodes of star formation occur mostly in ETGs in 
poor environments (e.g., [Thomas et al.|2010| ). 

While single-dish observations have the advantage of reach- 
ing good sensitivity in relatively short integrations, they lack the 
angular resolution needed to study the detailed H I morphology 
and kinematics. Starting from the 1980's an increasing number of 
galaxies was observed at higher angular resolution with interfer- 



ometers to perform such analysis (e.g., Shane 


1980|[Raimond et al. 


[1981| [Knapp & Raimond[[1984|[Sancisi et a 


l.[[ 1984 [van Gorkom 



et al.[ [l994 1995]). These and many later observations revealed 
a surprisingly large diversity of H I morphologies, ranging from 
very extended (tens of kpc), low column-density discs and rings 
to unsettled gas distributions indicative of recent gas accretion, gas 
stripping, or galaxy interaction and merging (see the review by [van[ 
[Gorkom & Schiminovich|T997l and [Hibbard et al.|200T] ). These re- 
sults revealed H I as a fundamental tracer of the assembly history of 
ETGs, and prompted work aimed at increasing the size of a sample 
with deep, homogeneous high-resolution H I data. 

A significant step in this direction was made by the combined 
study of ( [Morganti et al. 12006] herafter M06) and ( [Oosterloo et al.[ 
[20TO] herafter OlO) of 33 galaxies in the SAURON sample. Us- 
ing data taken with the Westerbork Synthesis Radio Telescope they 
detected about 10 percent and 2/3 of all galaxies inside and out- 
side the Virgo cluster, respectively (going down to H I masses of a 
few times 10^ M©). Exploiting the high resolution of their data they 
show that the H I is distributed on regular discs and rings in about 
half of the detections (1/3 of all galaxies), with radius ranging from 
~ 1 to tens of kpc (see also [Oosterloo et al.|2007[ ). 

Their analyses demonstrate that high sensitivity is crucial to 
detect the faintest signatures of H I accretion in ETGs, and that 
these are present in most detected galaxies. Furthermore, high an- 
gular resolution and the availability of multi- wavelength data make 
it possible to connect the H I to stars and multi-phase interstellar 
medium within the host galaxy. This is crucial to understand the 
role of H I in replenishing ETGs with cold gas, which could then 
fuel star- and disc formation. For example, M06 and OlO found that 
all ETGs surrounded by settled H I distributions host ionised gas 
within IRe, and that the two gas phases share the same kinematics. 
Furthermore, all galaxies with H I within \Re are detected in CO, 
and these systems are more likely to be detected in radio continuum 
too, indicating that some star formation is occurring (see OlO for a 
discussion of the origin of the radio emission). 

The main limitation of M06 and OlO is that the SAURON 
sample is small, hampering a statistically strong study of the re- 
lation between H I and other galaxy properties. For example, the 
number of galaxies was not sufficient to establish or rule out trends 
between H I and ETG dynamics or stellar populations, nor to study 
trends with environment beyond the simple Virgo vs. non- Virgo 
dichotomy. Also, the uncertainty on the fraction of ETGs hosting 
regular, rotating H I systems (undoubtedly a very interesting class 
of objects) was very large. 

Another limitation of the SAURON sample is that galaxies 
were selected to be evenly distributed on the MB-ellipticity plane 
rather than to follow the ETG luminosity function, so that the inter- 
pretation of the results in terms of galaxy evolution is not straight- 
forward. Works based on the HIPASS and ALFALFA surveys do 
not suffer from this limitation since both are blind surveys of large 
regions of the sky. However, they lack the resolution and sensitivity 
necessary to provide a picture as revealing as that emerging from 
M06 and OlO (and, in the case of ALFALFA, the sample studied 
so far is only 40% larger than the SAURON sample within their 
common Mb range). 

To overcome these limitations we extend here the study of 
M06 and OlO to a complete, volume-limited sample of 166 nearby 
ETGs (~ 5 times more galaxies) belonging to the ATLAS^^ sam- 
ple, while maintaining approximately the same sensitivity and res- 
olution of their observations. The large size and better selection of 
the sample studied here allow us to establish the detailed H I prop- 
erties of ETGs in the local Universe on a firm statistical basis. The 
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Figure 1. Distribution of Mk (left), Sio (middle) and (right) for galaxies in the ATLAS^^ sample (red line), ATLAS^^ H I sample (blue line) and H I 
SAURON sample (grey filled histogram). See the text for the definition of Eio and E3. 



availability of a large range of multi- wavelength data taken as part 
of the ATLAS^^ project is another element of continuity with work 
on the SAURON sample, and a major step forward relative to other 
studies. 

The aim of this article is to present the H I properties of galax- 
ies in the sample, how they depend on galaxy luminosity and en- 
vironment, and their relation to signatures of star formation in the 
host galaxy. Later work will explore how H I in ETGs relates to 
other galaxy properties derived from different datasets available 
within the ATLAS project, and the connection to simulations. 
This article is structured as follows. In Sec. [2] we introduce the 
galaxy sample. In Sec.|3]we describe radio observations and data 
reduction. In Sec.|4]we describe the H I morphology of the detected 
galaxies and introduce a classification scheme based on it. We also 
discuss the signatures of star formation in galaxies with different 
H I morphology, and the relation between H I and H2. In Sec. [5] we 
calculate upper limits on M(H I) of undetected galaxies. In Sec. [6] 
we discuss the Hi mass function of ETGs. In Sec. [7] we compare 
the distribution of M(H I), M(H i)/Lk and H I column density of 
ETGs and spirals. In Sec. [8] we investigate the relation between H I 
properties, galaxy luminosity and environment. Finally in Sec. [9] 
we summarise our findings and draw conclusions. 



2 SAMPLE 

We study the Hi properties of galaxies in the volume-limited 
ATLAS^^ sample, which includes 260 ETGs within 42 Mpc and 
brighter than Mk= -21.5 (Paper I). ETGs are selected from a par- 
ent sample as galaxies without spiral arms in available optical im- 
ages. 

Our Hi study is based on data taken with the Westerbork 
Synthesis Radio Telescope (WSRT). For observability reasons we 
select only the 170 galaxies with declination 6 > \0 deg. We 
also exclude all 4 galaxies closer than 15 arcmin to the centre of 
the Virgo galaxy cluster (NGC 4476, NGC 4478, NGC 4486 and 
NGC 4486A; observations close to Virgo A do not provide H I data 
of sufficient quality). Therefore, we study here a sub-sample of 166 
ATLAS^^ ETGs, 61 percent of the total. We refer to this as the 
ATLAS^^ H I sample. Of these galaxies, 39 reside inside the Virgo 
galaxy cluster and 127 outside it (Paper I). All 166 galaxies are 



listed in Table |B1| together with properties derived in the present 
work. 

Figure [1] shows the distribution of galaxy absolute magnitude 
Mk, large-scale environment density Eio and local environment 
density E3 for the ATLAS^^ H I sample and the full ATLAS^^ sam- 
ple. These parameters are given in Paper I (M^) and in Paper VII 
(Sio and S3; these are defined as the number density of galaxies 
contained within a 600-km/s-deep cylinder whose radius is equal to 
the distance from the tenth and third closest galaxy, respectively). 
Fig. [1] shows that the distribution of these parameters is essentially 
the same for the ATLAS^^ H I sample and the full ATLAS^^ sam- 
ple. 

This H I survey expands the study performed by M06 and OlO 
on the SAURON sample. All but one of the 33 galaxies studied 
by M06 and 010 belong to the ATLAS^^ sample. Figure [l] shows 
the distribution of galaxies common to the two samples. SAURON 
galaxies are not distributed very differently if one takes into account 
statistical uncertainties. The main improvement of the present study 
are the selection of the sample (which, as discussed in the previous 
section, is unbiased in galaxy luminosity) and sample size, so that 
we now have a more representative sampling of the ETG properties. 



3 Hi DATA 

In this section we describe radio observations and data reduction. 
For some galaxies data were taken as part of earlier studies. Refer- 
ence to the corresponding papers is given in the text below and in 
TablelBTI 

3.1 Observations 

We follow two different observational strategies for galaxies in- 
side and outside the Virgo cluster. We observe all but one of the 
127 galaxies outside Virgo with the WSRT (the only exception is 
NGC 4762, for which we use ALFALFA data). WSRT data were 
taken as part of previous projects for 18 of these 127 objects ( |J6zsa| 
let al.|2009| M06 and OlO). Another 3 fall within the field of view of 
these earlier observations. All remaining 105 galaxies are observed 
for ATLAS^^. Both old and new data are taken using the WSRT in 
the maxi-short configuration. We observe using a bandwidth of 20 
MHz sampled with 1024 channels, corresponding to a recessional 
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velocity range of ~ 4000 km/s and a channel width of ~ 4 km/s 
for the H I emission line. The only exception is NGC 2685, ob- 
served by |J6zsa et aL| ( |2009] ) with 1024 channels over a 10-MHz 
bandwidth. 

We observe each galaxy for 12 h. For some galaxies, earlier 
data are deeper than our 12 h integration (see Table |B1| column 8) 
and we use those in our analysis. 

Data are less homogeneous for the 39 galaxies inside Virgo. 
We take WSRT data with the same array configuration and cor- 
relator set-up described above for 3 galaxies. Equivalent data are 
available for 13 more galaxies in OlO. A field including 2 more 
ATLAS^^ galaxies was observed with the WSRT by Oosterloo and 
collaborators with the same WSRT configuration and correlator set- 
up described above, and we use their data in our analysis. We use 
the ALFALFA database to look for H I emission in the remaining 
21 Virgo galaxieQ We find 2 galaxies with a possibly associated 
ALFALFA H I detection: NGC 4694 and NGC 4710. We make use 
of VLA data taken for the VIVA project to study the former ( Chung 
|et a l. 2009 ); and we observe the latter with the WSRT using the 
same array configuration and correlator set-up as above. We use 
ALFALFA spectra to derive M(H I) upper limits for the remaining 
19 galaxies (see Sec.|5]). 

3.2 WSRT data reduction and products 

We reduce the WSRT data in a standard way using a dedicated 
pipeline based on the Miriad reduction package ( |Sault, Teuben| 
|& Wright] 1 1995| ). The standard pipeline products are Hi cubes 
built using robust=0 weighting and 30-arcsec-FWHM tapering (see 
[Brig gs 1995 for an explanation of the robust parameter). The dis- 
cussion below is based on the analysis of these cubes unless other- 
wise stated. The pipeline is validated by comparing its products to 
those obtained with manual reduction by M06 and OlO for galaxies 
in common. 

The typical angular resolution of the cubes is 35 x 35/ sin (S) 
arcsec^, where 6 is the declination. All cubes are made at a velocity 
resolution of 16 km/s (after Manning smoothing). The noise ranges 
between 0.4 and 0.7 mJy/beam. The median noise is 0.5 mJy/beam, 
and 90 percent of the cubes have noise below 0.6 mJy/beam. At 
the median noise level the 5<t column density threshold within one 
velocity resolution element is ~ 3.6 x 10^^ x sin ((5) cm~^. The 
ATLAS^^ H I sample covers the 6 range 10-60 deg, corresponding 
to A^(H I) in the range 0.6 - 3.1 xlO^^ cm'l 

For each galaxy we build a total H I image by summing along 
the velocity axis all emission included in a mask constructed as 
follows. A pixel belongs to the mask if the absolute value of its 
flux density is above 4cr in at least one of the following cubes: (/) 
the original cube; (ii) the three cubes obtained smoothing the orig- 
inal cube in velocity with a Hanning filter of FWHM 32, 64 and 
112 km/s; ( the original cube smoothed with a 60-arcsec-FWHM 
Gaussian beam; (iv) the three cubes obtained by Hanning smooth- 
ing the 60-arcsec cube as in point ( ii). This allows us to be sensitive 
to H I emission over a wide range of angular and velocity scales. 
The binary masks are then enlarged by convolution with the synthe- 
sised beam, and only mask pixels whose value is above 50 percent 
of the convolved mask peak are kept. This removes a large fraction 
of spurious detections with small angular size. It also ensures that 
we do not miss faint gas just below the detection threshold at the 
edge of the H I distribution. 



http : //arecibo . tc . Cornell . edu/ hi archive/ alfalfa 



Total H I images are obtained as the zero-th moment of the 
masks built following the above procedure. These images contain 
also pixels with negative surface brightness because we apply the 
detection algorithm to the absolute value of the flux density. We 
use these "negative" pixels to estimate the column density detec- 
tion threshold of the images. In practice, we determine the surface 
brightness threshold -5 th below which only 5 percent of the neg- 
ative pixels are retained. The final total H I image includes only 
pixels with \S\ > +»Sth. The typical value of Sth corresponds to an 
H I column density of a few times 10^^ cm~^, the exact value vary- 
ing from image to image. 

We also build H I cubes with robust=l weighting (close to nat- 
ural weighting) and no tapering. Below we refer to this weighting 
scheme as i?01. The noise in R0\ cubes is a factor of ~ 1.5 lower 
than in the standard cubes. However, its pattern is very patchy and 
the overall image quality is lower. Therefore, we use these cubes 
only when the detection algorithm described above reveals low- 
surface-brightness emission missing from the standard cubes (see 
Table |Bl| column 7). In these cases, the H I mass given in Table |Bl| 
is derived from the R0\ cube. 

We note that the ROl cubes have slightly better angular reso- 
lution than the standard cubes: 25 x 25 /sin(6) arcsec^. Therefore, 
the formal column density sensitivity is ~ 30 percent worse. 



4 H I MORPHOLOGY AND ITS RELATION TO THE 
HOST GALAXY 

We detect 53/166 = 32 percent of all ETGs in the ATLAS^^ Hi 
sample. Unlike the detection rate presented in M06 and OlO, this 
does not include a number of cases where H I is detected only in 
absorption or where we cannot securely assign H I to the observed 
ETG because of confusion with neighbouring galaxies (see below). 

The H I detection rate depends strongly on environment den- 
sity, being 4/39=10 percent inside Virgo and 49/127=39 percent 
outside it (see Sec. |8] for a comparison to previous results). We 
postpone a full discussion of environmental effects to Sec. |8] In 
this section we describe the H I morphology of all detected galax- 
ies regardless of their environment. 

4.1 H I morphological classes 

We show total H I images of all detected ETGs in Fig. |Al| The fig- 
ure suggests the existence of a few broad types of H I morphology, 
which we describe with the following H I morphological classes: 

• D (large discs) - most of the H I rotates regularly and is dis- 
tributed on a disc or ring larger than the stellar body of the galaxy. 

• d (small discs) - most of the H I rotates regularly and is dis- 
tributed on a small disc confined within the stellar body of the 
galaxy. 

• u (unsettled) - most of the H I exhibits unsettled morphology 
(e.g., tails or streams of gas) and kinematics. 

• c (clouds) - the H I is found in small, scattered clouds around 
the host galaxy. 

After examining morphology and kinematics of the detected H I we 
assign each galaxy to one of these classes. The H I class of detected 
objects is given in Table |Bl| and indicated in Fig. |Al| Further notes 
on the H I morphology and kinematics of individual galaxies are 
given in Table |Bl| column 7. 

This classification scheme is an attempt to simplify the vari- 
ety of H I morphologies seen in Fig. |A1| It also reflects our view 
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Table 1. Number of galaxies in the various H I classes as a function of Virgo 
membership (see Sec.jS} and the presence of CO and dust/blue features (see 
Sec.|44). 



Hi 


all 


inside 


outside 


CO 


dust/blue 


class 


galaxies 


Virgo 


Virgo 


det. 


features 


all 


166 


39 


127 


38 


36 


undet. 


113 


35 


78 


15 


13 


det. 


53 


4 


49 


23 


23 


D 


24 


1 


23 


10 


10 




(10) 


(1) 


(9) 


(1) 


(2) 


d 


10 


1 


9 


8 


10 


u 


14 


2 


12 


5 


3 




(4) 


(0) 


(4) 


(0) 


(0) 


c 


5 





5 









Numbers in parenthesis indicate the number of galaxies within that class 
for which H I is not found in the centre of the host galaxy. All J's and none 
of the c's host H I in their centre. 

that different H I morphologies give different indications on the as- 
sembly and gas-accretion history of the host galaxy. For example, 
D's and J's are galaxies with a relatively quiet recent history on a 
time- scale proportional to the orbital time of the rotating gas (this 
is obviously larger in D than d systems; as an example, it is ~ 1 Gyr 
for the D galaxy NGC 6798, and ~ 200 Myr for the d NGC 5422). 
On the contrary, galaxies in class u, where most of the gas is not ro- 
tating regularly around the stellar body, have recently experienced 
tidal forces or episodes of gas accretion/stripping involving a large 
fraction of their gas reservoir. In other words, their H I content is 
evolving, and for some of them the current H I properties may be 
just a transient phase. In contrast, in D objects the H I has likely 
been part of the galaxy (and orbiting around the stellar body) for 
many gigayears. Finally, galaxies in class c are surrounded by H I 
not obviously associated to them, and we regard them as objects 
with similar H I properties as undetected ETGs. 

Similarly, the distinction between D and d galaxies, which is 
entirely based on the size of the H I distribution, is motivated by our 
view that they are fundamentally different objects (see Sec. |4.3| ). In 
order to make this classification as objective as possible we define 
the H I radius /?(H I) as the maximum distance of the N{Y{ l)=5 X 
10^^ cm~^ isophote from the galaxy centre (deconvolved with the 
beam size). We then compare i?(H I) to the optical effective radius 
Re given in Paper I. We adopt i?(H l)= 3.5 x as the dividing line 
between small and large H I systems (for a de Vaucouleur profile 
this radius includes 90 percent of the total light) 

Given this classification scheme, D galaxies have M(H I) in 
the range 10^-10^^ M©, and the gas extends out to many tens of 
kpc from the centre of the galaxy (in some cases more than 10 
Rq). In contrast, M(Hl) in d systems is typically below 10^ M©, 
and the size of the gas distribution is up to a few kpc. Galaxies in 
group u contain between a few times 10^ and 10^^ M© of H I, often 
stretching to many tens of kpc from the galaxy. Finally, galaxies in 



^ The only galaxy whose classification would not fit in our picture is 
NGC 3414, which we regard as hosting a large Hi system although its 
7?(H i) is smaller than 3.5x7?^. However, the gas in this galaxy is distributed 
on a polar configuration, so that it may be more appropriate to compare 
/?(Hi) to X (1 - e) where e is the optical ellipticity given in Paper 11. 
When this is done the ratio between H I and optical radius is larger than 3.5 
so that we classify this galaxy as D. 



group c have M(H I) below a few times 10^ M©, comparable to the 
H I mass upper limits described in Sec. [5] 

4.2 A continuum of H I morphologies 

We find 24, 10, 14, and 5 galaxies in group D, d, u and c, respec- 
tively, corresponding to 14, 6, 8 and 3 percent of the sample. Only 
I D, \ d and 2 w's belong to the Virgo cluster. Outside Virgo 18, 
7, 9 and 4 percent of all galaxies belong to group D, d, u and c, 
respectively. About 40 percent (10/24) of all D's are rings (or discs 
with a central H I hole). We summarise these results in Table [T] 

The above figures show that within the ATLAS^^ volume as 
many as 1/5 of all ETGs (and 1/4 of all ETGs outside Virgo) host 
rotating H I distributions (groups D and d). These objects represent 
the majority of our detections, 64 percent, in agreement with M06 
and OlO. Therefore, ifH I is detected in an ETG it will most likely 
be distributed on a disc or a ring. 

Within the adopted classification scheme it is not always easy 
to assign a galaxy to a given class. In fact, we find a number of 
objects intermediate (or in transition) between classes. For the sake 
of simplicity we will make use of this classification in the rest of 
the article, but the ATLAS^^ H I sample reveals rather a contin- 
uum of H I morphologies going from settled, rotating systems to 
progressively more and more disturbed ones, and from the latter to 
systems of scattered H I clouds which may or may not be related 
to the nearby ETG. We guide the reader through this continuum of 
morphologies in the rest of this section. 

Figure[2] shows a sequence of increasingly less regular H I dis- 
tributions (left to right) whose size is large relative to that of the 
host galaxy. The sequence shows the continuity of H I morphology 
going from the most regular D objects (left) to the most irregular, 
large w's (right). It begins on the left with an H I ring, likely a res- 
onance related to a bar (see also the stellar ring coincident with 
the H I distribution). The morphology (top row) and kinematics 
(bottom row) of the H I ring are extremely regular. The sequence 
continues on the second panel with a disc whose rotation and mor- 
phology are somewhat less regular. In this galaxy a tail of gas is 
visible to the west, and a small neighbour to the east might be in- 
teracting with our target. Although the H I kinematics exhibits an 
overall ordered rotation (but misaligned relative to the stellar kine- 
matics), it is clear that some of the gas has not settled yet within the 
galaxy potential (or has been disturbed recently). Systems like this 
provide a link between Z)'s and galaxies in group w, where most 
(but not necessarily all) of the gas is unsettled. 

The third panel shows a galaxy that we classify as u. However, 
there are hints that at least part of the H I is rotating around the 
stellar body (see the position-velocity diagram on the bottom row 
and the description of this galaxy in |Sancisi et al.|1984] and M06). 
Systems like this could move towards the D class if their evolu- 
tion proceeds without significant merger/accretion events for a few 
gas orbital periods. Finally, the fourth panel shows completely un- 
settled H I ( jPuc et al.|2011| herafter Paper IX, discuss the optical 
signature of a recent merger in this galaxy). It is possible to see 
the sequence in Fig. [2] as one of decreasing time passed since the 
last major episode of gas accretion or stripping (relative to the gas 
orbital period). 

A similar sequence can be drawn for H l-detected galaxies 
where the size of the gas distribution is comparable to that of the 
stellar body. Such a sequence, shown in Fig. [3] illustrates the con- 
tinuity of H I morphology going from the most regular d objects 
(left) to the most irregular w's (right). The sequence starts with a 
very regular small H I disc. It continues with a very faint, small 
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Figure 2. A sequence of H i-rich ETGs with increasingly less regular gas configurations (left to right). The sequence shows galaxies with very large H I 
distributions relative to the stellar body. The top row shows total-H I contour images. We refer to the caption of Fig. |Al| for a description of the content of each 
image. In this figure the top-right scale bar indicates 10 kpc at the galaxy distance. The bottom row shows position- velocity diagrams of galaxies in the top 
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drawn at 1.0 x 2" mJy beam~^ ,n - 0, 1, 2, .... The cross on the bottom left indicates 10 kpc along the horizontal axis and 50 km s~^ along the vertical axis. 



disc connected spatially and in velocity to a cloud outside the stel- 
lar body. This system appears fairly regular within the stellar body 
(given the low signal-to-noise ratio), but the outer cloud indicates 
that some accretion may be on-going. In fact, the H I cloud sug- 
gests that small H I discs may form by accreting gas from small 
companions or the surrounding medium (OlO). 

The third panel in Fig. [3] shows a system which is difficult 
to classify. We detect H I only on one side of the galaxy, and be- 
cause of this we classify it as u. However, the position-velocity di- 
agram shows a hint that gas may be present also on the other side 
of the galaxy at opposite velocity relative to systemic. This is the 
behaviour expected for a disc, so this galaxy may actually be a mis- 
classified d. We conclude the sequence with a system where H I is 
detected far from the stellar body but is clearly connected to it. The 
modest amount of gas and the size of the H I distribution relative to 
the stellar body suggests that this system may evolve towards the d 
class if left undisturbed. 

Another step in this continuum of H I morphologies is illus- 
trated in Fig.|4] In this figure we show that some the unsettled sys- 
tems appear as a bridge between class u and class c. The sequence 
starts from an extended, unsettled H I system. The large number 
of gas-rich neighbours suggests that galaxy interaction may be re- 
sponsible of the disturbed H I morphology (see also [Verheijen &| 
|Zwaan|2001| ). The second and third galaxies on the sequence are 
surrounded by a few H I clouds forming a coherent gaseous system 
- in both cases the clouds are connected to each other and to the 
ETG in velocity and are likely to be the densest clumps of an H I 
tail (for NGC 5557 see Paper IX). These systems represent a natu- 
ral link to galaxies in group c (represented by the fourth object in 



the sequence), where scattered clouds are detected around an ETG 
but their connection to it is much less obvious. 

The final piece in this continuum of morphologies is provided 
by the continuity between c and undetected objects. Many c galax- 
ies are difficult to classify as it is sometimes not clear whether we 
are seeing the densest gas clumps of a w system, or gas is instead 
floating unbound in the inter-galactic medium. In fact, a number 
of ETGs which we classify as non detections are surrounded by H I 
with properties similar to those of c galaxies. We show these objects 
in Fig. |A2| (the figure includes also galaxies close to or interacting 
with an H l-rich neighbour, so that confusion prevents us from es- 
tablishing whether they contain any H I). This ambiguity has little 
effect on our discussion in the following sections as the H I mass 
of c objects is comparable to the detection limit of our observations 
(see Sec.|5]). 



4.3 Small and large H I discs: relation to the kinematics of 
the host galaxy 

At the beginning of this section we introduced a criterion to distin- 
guish between large and small rotating H I systems (i.e., D and d). 
We consider this distinction important as there are fundamental dif- 
ferences between these objects (besides their size). As mentioned, 
D's contain between 10^ and M© of H I, while in J's M(H I) is 
typically below 10^ M©. Furthermore, the H I kinematics is aligned 
with the stellar kinematics in most galaxies belonging to group d - 
8/10 objects. Figure |5] shows position- velocity diagrams of all J's 
along the stellar kinematical major axis given in Paper II (perpen- 
dicular to it in the case of NGC 3499 - see below). Comparing 
these diagrams with the stellar velocity fields published in Paper II 
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we see that the H I rotation is aligned with that of the stars in all 
galaxies except NGC 3032, where the H I is counter-rotating, and 
NGC 3499, where the H I is polar. Even in these two cases a con- 
nection between H I and stellar body exists. In the former the gas 
disc is aligned with a kinematically-decoupled stellar component 
( |McDermid et al.|2006| OlO). In the latter the H I is aligned with a 
faint stellar disc and a dust lane visible in the SDSS image. 

In contrast with this result, visual inspection of the stellar ve- 
locity fields (Paper II) and the HI cubes shows that misaligned 
discs/rings are found in more than half (14/24) of all D galax- 
ies. This includes a number of polar rings and counter-rotating 
discs/rings (see Fig. |Al| and column 7 in Table |B1| ). From a kine- 
matical point of view, it seems that d H I distributions are tightly 
coupled to their host galaxy, while in D's the H I kinematics is of- 
ten different from the stellar one. A future paper will address the 
relation between H I and stellar kinematics in more details, includ- 



ing the classification of ETGs in the ATLAS^^ sample in slow- and 
fast rotators (Paper III). 



4.4 H I, H2 and star formation signatures 

As part of their study of 33 ETGs belonging to the SAURON sam- 
ple, 010 show that all galaxies where H I is detected within ~ IR^ 
contain molecular gas in their centre, and that H2 dominates the 
ISM in this region. On the contrary, the CO detection rate is ~ 20 
percent for ETGs with no central H I detection. Here we revisit their 
result by comparing the H I properties of ETGs in the ATLAS^^ H I 
sample to results of Paper II (which provides a catalogue of dust 
discs, dusty filaments and blue regions in all ETGs in the ATLAS 
sample) and Paper IV (where we discuss the ATLAS^^ CO survey). 

We detect dust/blue features and CO in 36 and 38 of the 
166 ETGs in the ATLAS^^ Hi sample, respectively. The pres- 
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Figure 5. Position- velocity diagrams of all d galaxies along the stellar kinematical major axis given in Paper II (perpendicular to it for NGC 3499). In all 
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ence of dust/blue features and molecular gas are very tightly re- 
lated to one another (Paper IV), so that most galaxies with CO also 
have dust/blue features and vice- versa. With reference to Table [T] 
dust/blue features and CO are present in 43 + 9 percent of all H I 
detections. In contrast the dust/blue and CO detection rate is 13 + 3 
percent for the H I non-detections. Therefore, galaxies detected in 
H I are ~ 3 times more likely to contain signatures of star formation 
than those with no detectable neutral hydrogen. 

We make use of the classification presented above to investi- 
gate how H I morphology is related to the occurrence of star for- 
mation. In the previous section we discuss the kinematical link 
between H I and stellar body in d ETGs. This result is comple- 
mented by the detection of star formation signatures in all these 
galaxies. As illustrated in Table [T] dust/blue features are detected 
in 10/10 = 100 + 32 percent of thenj^ Furthermore, we detect CO 
in 8/10 = 80 + 28 percent of them. The exceptions are NGC 3499 
and NGC 5422, for which the upper limit on M(H2) is ~ 5 x 10^ 
Mo. 

A further indication that stars are being formed efficiently in 
d galaxies comes from the analysis of the ISM composition in their 
central region. As done in OlO, we compare the H2 mass given 
in Paper IV (measured within the ~20 arcsec IRAM beam) to the 
H I mass contained in one WSRT beam (whose minor axis is ~ 30 
arcsec) centred on the galaxy. This comparison shows that in J's the 
M(H2)/M(H I) ratio is high - always above ~ 1, and > 10 in 5 out 
of 8 galaxies detected in CO. Therefore, the ISM in these galaxies 
has similar properties as that in the central region of spirals ( |Leroy| 
|et al.|200"8| ), confirming that in many of them H I is being turned 
into H2 efficiently. 

Signatures of star formation are less common in D galaxies 



^ Note that NGC 4150 is not Hsted as containing dust/blue features in Paper 
II because of the poor SDSS image quality. However HST observations by 
[Crockett et al.lpOll) show that also this galaxy is characterised by a dusty 
and star-forming core. 



and depend on whether the H I is distributed on a ring (10 objects 

- see Table [T]) or extends all the way to the centre of the galaxy (14 
objects). We detect CO in 9 of the ETGs with central H I (64 + 21 
percent) and dust/blue features in 8 of them (57 + 20 percent). In 
contrast only 1 galaxy with an Hi ring is detected in CO (10 + 
10 percent) and 2 (including the CO detection) exhibit dust/blue 
features (20 +14 percent). 

The physical state of the ISM in these objects may be differ- 
ent from that in J's as the M(H2)/M(H I) ratio is somewhat lower. 
When CO is detected, we find M(H2)/M(H I) values in the range 0.5 

- 3 for 80 percent of the Z)'s, and never larger than ~ 10. When CO 
is not detected the upper limit on M(H2)/M(H I) is always above 
~ 1, so that these systems are not necessarily different from those 
with a CO detection (they may simply contain too little central H I 
for molecular gas to be detected). Therefore, also in D's most of 
the central ISM is found in the form of molecular gas, but the con- 
version of H I to H2 may be less efficient than in d's. One possibile 
explanation is that H I in the centre of these galaxies has lower col- 
umn density. 

A similar result holds for u galaxies. These systems are un- 
settled so that it is not always easy to judge whether H I is present 
within the stellar body, or is simply projected onto it on the sky. 
However, Hi is certainly not present in the central region of 
NGC 3193, NGC 4026 (a ring with a large and massive H I tail), 
NGC 5198 and NGC 5557. None of these systems hosts CO or 
dust/blue features. Of the remaining 10 u galaxies, 3 host dust/blue 
features (30 +17 percent) and 5 host molecular gas (50 + 22 per- 
cent). As in D's, the typical (upper limit on) M(H2)/M(H I) is of the 
order of unity, and none of these objects has M(H2)/M(H l)> 10. 

Consistent with the above discussion, none of the 5 c's exhibits 
dust/blue features or is detected in CO. We therefore confirm OlO 
result that the detection rate of star-formation signatures is high 
(~ 2/3) for ETGs hosting neutral hydrogen within ~ \Rq, so that 
H I seems to provide material for star formation. In particular, 100 
percent of all J's, ~ 60 percent of all D's with central H I, and ~ 50 
percent of all w's with central H I exhibit such signatures. At the 
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Figure 6. Peak signal to noise (S/N)inax of all galaxies detected in the standard, 12-h H i cubes (see text) plotted against the size of the H i system in units of 
the beam (left), and the spectral line width at the position of (S/N)jnax (right). The top H symbols represent the FWHM of the Manning filters used to smooth 
the H I cubes when searching for emission (see Sec.jSj. 



sensitivity of our data, even galaxies with central H I but no star 
formation signatures are consistent with the picture that the central 
ISM of ETGs is dominated by molecular gas - they may simply 
contain too little H I for CO to be detected. In contrast with these 
results, galaxies surrounded by an H I distribution with a central 
hole, and galaxies with no detectable H I, show signs of recent star 
formation in a minority of the cases (~ 15 percent). 

These trends indicate that H I may play an important role in 
fuelling rejuvenation and the growth of a stellar disc in ETGs. In 
particular, in a number of J's (the ETGs with the highest H l-to-H2 
conversion rate in their centre) the H I morphology suggests that 
recent accretion (possibly from small companions) is the source of 
the ISM (e.g., NGC 3032, NGC 3489, NGC 4150 - see OlO). Yet, 
15/38 CO detections and 13/36 ETGs with dust blue/features are 
not detected in H l. The H2 mass of these systems is not particu- 
larly low, so they are galaxies with a very high M(H2)/M(H I) ratio 
(lower limits on this are above ~ 5 in most of these cases). Hosting 
H I at the present time (and at the sensitivity of our observations) is 
therefore not a necessary condition for an ETG to host star forma- 
tion. It is possible that some ETGs have been stripped of most their 
H I but not their molecular gas if they live in a dense environment 
with a hot medium (see Sec. [8] and Paper IV). However, more than 
half of the galaxies with signs of star formation but no detected H I 
live outside Virgo, so this cannot be the only explanation. These 
results highlight the complex gas accretion history and interplay 
between different ISM phases of ETGs. 



[Sadler, Jenkins & Kotanyi|1989||Mauch & Sadler|2007] ). A com- 
plete discussion of these results will be presented in a future paper. 

We look for H I absorption in all objects with a continuum 
flux above 8 mJy. This results in three detections of objects where 
H I is observed only in absorption (NGC 5322, NGC 5353 and 
PGC 029321), while in three more objects absorption is detected 
in addition to emission (NGC 2824, NGC 3998 and NGC 5866; 
see Table |Blj ). The continuum flux density of objects detected in 
H I absorption ranges between 8 and ~ 30 mJy. The H I absorption 
is relatively narrow, 50 - 80 km/s, and in all cases it is centred on the 
systemic velocity. The absorbed flux ranges between 3 and 4 mJy. 
We derive optical depths of 7 - 10 percent and column densities of 
4-1 X 10^0 cm-2 (assuming T^pin = 100 K). 

This absorption may be the result of small disc structures (i.e. 
discs that are not large enough to be seen in emission if a relatively 
strong continuum source is present). For example, two of the three 
galaxies detected only in absorption exhibit dust or blue features 
within the stellar body, consistent with the result described above 
for d galaxies. Only one of them is detected in CO (PGC 029321). 

It is worth noting that only 19 objects among the ~ 140 ob- 
served with the WSRT have radio continuum flux density S \aghz > 
8 mJy. The detection rate of H I in absorption is therefore ~ 30 
percent for galaxies with sufficient continuum emission. This un- 
derlines that the H I detection rate of our study (in emission and 
absorption) is just a lower limit to the fraction of ETGs containing 
neutral hydrogen gas. 



4.5 H I in absorption 

We conclude this section by noting that H I gas can also be detected 
in absorption against a radio continuum central source (1.4 GHz 
continuum images are produced as part of the pipeline described in 
Sec.[3]). 

For each object, we obtain continuum images using the line- 
free channels. The typical noise of the images is ~ 60 juJy beam~^ 
and more than 40 percent of the galaxies were detected with ra- 
dio power in the range 10^^-10^^ W Hz~\ consistent with previ- 
ous studies of the continuum emission of this kind of objects (e.g.. 



5 M(H I) DETECTION LIMITS 

The analysis presented in the following sections is based on the 
images in Fig. |A1| the classification discussed above, the H I mass 
measured from these images and listed in Table |B1| as well as on 
M(H I) upper limits for the undetected galaxies. In this section we 
discuss how we calculate these upper limits. 

The data cubes described in Sec.[3]give the H I surface bright- 
ness distribution projected on the sky per unit velocity interval 
along the line of sight. The detection of H I in a galaxy is there- 
fore driven by its peak surface brightness relative to the noise level 
at the angular and velocity resolution of the data, and not by the 
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Figure 7. Distribution of the M(H i)/M(H ratio for galaxies detected 
in the standard, 12-h cubes (see text). 

total mass of H I hosted by the galaxy. This makes the calculation 
of upper limits M(H l)iim for undetected galaxies rather delicate. 

The typical approach in the literature is to calculate M(H l)iiin 
as the minimum detectable H I mass within one angular-resolution 
element centred on the galaxy, and integrating over a velocity in- 
terval representative of the gas line-width. For example, if we 
calculate M(H l)iiin as a 3 cr signal within one synthesised beam 
and within a 200 km/s velocity interval, we obtain M(Hl)iim= 
2.0 • (J/10Mpc)2 X 10^ Mq. However, a larger Hi mass can be 
hidden below the noise if it is spread over many beams. The first 
question we wish to answer is therefore whether we are missing a 
population of extended, massive H I distributions with column den- 
sity below the noise of our observations. 

To find an answer we calculate the maximum signal-to-noise 
ratio (S/N)inax of the H I emission for galaxies detected in the stan- 
dard cubeQ This calculation is made including data-cubes at all 
angular and velocity resolutions listed in Sec 3 (i.e., all possible 
combinations of the velocity resolutions 16, 32, 64 and 112 km/s, 
and of the angular resolutions ~ 30 and 60 arcsec). The value of 
(S/N)inax indicates how easily a galaxy is detected above the noise. 

The left panel in Fig. [6] shows (S/N)inax plotted versus the 
area A(H I) occupied by the detected gas in units of the beam area 
A(beam). This plot shows a trend of increasing (S/N)inax with in- 
creasing A(H l)/A(beam). In particular, the majority of the detec- 
tions have very high (S/N)jnax (above 10), and this is true for nearly 
all detections larger than 10 beams. There is no hint of a population 
of very extended H I systems close to the noise of our observations. 

The second question we wish to answer is what angular size 
and line- width we should assume to calculate M(H l)iiin for the un- 
detected galaxies. To find an answer we study the distribution of 
A(H l)/A(beam) and of the line-width W of low-(S/N)inax galaxies 
((S/N)jnax < 10)- The typical low-(S/N)jnax detection is rather small 
on the sky. For these systems, the mode of the A(H l)/A(beam) dis- 
tribution is 1, and 2/3 of them have an area smaller than 6 beams, 
which we adopt as angular size for the calculation of upper limits. 

^ We therefore leave out NGC 3499, NGC 4150 and NGC 7332, which are 
detected only in the ROl cubes, and NGC 4406, for which only a total H I 
image (and not the cube, necessary for this analysis) is available. We note 
that for some galaxies the available cubes are obtained after multiple 12-h 
integrations (see Table [Bl} . For consistency, we re-analyse a single 12-h 
integration for these galaxies, and use the result in this section. 



Figure 8. Normalised histogram of M(H i) for ETGs in the ATLAS^° par- 
ent sample. The dashed red line refers to H i-detected ETGs. The solid red 
line is the corrected histogram (see text). Error bars are obtained assuming 
a binomial distribution. The shaded histogram represents the distribution 
of upper limits for undetected galaxies. The solid black line represents the 
best-fitting Schechter function. The dotted black line represents a Schechter 
function which is equal to the best-fitting one except for the slope a, fixed 
to the value obtained by Zwaan et^p005a) - 



We define W as the line width including 80 percent of the 
detected emission at the sky coordinates of (S/N)^^^^. The right 
panel in Fig. [6] shows that (S/N)jnax has little relation to W. The 
top H symbols represent the FWHM of the Hannning filters used 
to smooth the H I cubes in the spectral direction. They demonstrate 
that we have properly sampled the W range covered by the data. 
We find that median and mode of the overall W distribution are 
~ 50 km/s. Three quarters of all low-(S/N)jnax galaxies have W be- 
low this value, which we adopt as line- width for the calculation of 
upper limits. 

Based on these results, we define M(H l)iiin as a 3<t signal ob- 
tained integrating over 6 beams and 50 km/s. For the typical noise 
of 0.5 mJy/beam at a 16 km/s velocity resolution, this corresponds 
to M(Hl)iiin= 2.4 • (J/10Mpc)2 X 10^ M©. To verify that this is 
a reasonable upper limit we show in Fig. |7] the histogram of the 
M(H l)/M(H l)iiin ratio for the detected galaxies. For consistency, 
we make use of the M(H I) value measured from the standard, 12 
h-integration cube for all galaxies (this is lower than the M(H I) 
value reported in Table |B1| for galaxies were additional Hi emis- 
sion is detected in ROl or deep cubes). Only 6 of the 49 galaxies 
used for this analysis have M(H I) below M(H 1)11^ (NGC 2824, 
NGC 3182, NGC 3608, NGC 4710, NGC 5422, NGC 5866). This 
confirms the validity of our definition of upper limits. 

As mentioned in Sec. [3] we use ALFALFA spectra to deter- 
mine M(H l)iim for about half of all Virgo ETGs. Upper limits from 
ALFALFA data are calculated on one resolution element. The AL- 
FALFA beam has a FWHM of ~ 3.5 arcmin. Therefore, its area 
is larger than the 6 WSRT beams over which we calculate upper 
limits for the WSRT data. As above, we assume a line-width of 
50 km/s. ALFALFA spectra have a typical noise of ~ 3 mJy/beam 
at 11 km/s resolution, so that M(Hl)iim= 5.0 • (d/lO Mpcf X 10^ 
Mq. The noise value of individual spectra was kindly provided by 
Riccardo Giovanelli. 
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6 H I MASS FUNCTION 

The large size of the ATLAS^^ H I sample and the good sensitiv- 
ity of our WSRT observations result in a relatively large number 
of H I detections, making it possible to study the H I mass function 
of ETGs. Figure |8] shows the distribution of M(H I) for ETGs in 
the ATLAS^^ H I sample. The red dashed line shows the observed 
distribution of H l-detected ETGs. The shaded histogram shows the 
distribution of upper limits for undetected galaxies. Both distribu- 
tions are normalised by the total number of galaxies in the sample. 

We correct the histogram of detections for the incompleteness 
of our observations at low M(H I). We do this by normalising the 
number of galaxies detected in each bin by the number of galaxies 
detectable in that bin (i.e., the number of galaxies with M(H l)iim 
below the upper end of the bin). This is similar to the V/Vmax cor- 
rection typically done in large flux-limited surveys ( [Schmidt|1968| ). 
We show the corrected histogram with a red solid line in the figure. 
Error bars in the figure assume a binomial distribution. 

The correction can be applied only to a sample of galaxies 
each with M(H l)^M(H 1)11^. Three H l-detected ETGs do not meet 
this criterion (NGC 3182, NGC 3499 and NGC 7332). The cor- 
rected histogram in Fig. [8] are drawn ignoring these galaxies (i.e., 
considering them effectively as undetected). 

The corrected M(H I) distribution of ETGs is flat between a 
few times 10^ and ~ 3x10^ M©, and declines rapidly above the lat- 
ter value. At lower masses it exhibits a slight decline all the way to 
M(H l)~ 10^ Mq, where it may start rising again (but uncertainties 
in this mass range are large). We quantify this behaviour by fitting a 
Schechter function to the corrected M(H I) distribution ( |Schechter| 
[1976): 

dn oc (M/M*y^'' Jlogio (1) 

where dn is the number of galaxies within a logarithmic bin of 
width J log 10 M. The solid black line in the figure represents the 
best fit, which has parameters M* = (1.8 + 0.7) x 10^ M© and 
a = -0.68 + 0.16. 

The above mass function is derived neglecting Hi non- 
detections. Figure [8] shows that this may lead us to underestimate 
the mass function (and overestimate a) below M(Hl) ~ 5 x 10^ 
Mq. We therefore repeat the fit using only the distribution above 
this mass. The result is a = -0.79 + 0.23, consistent with the fit on 
the full M(H I) range. We note that the shape of the ETG H I mass 
function is not driven by environmental effects. The same shape is 
recovered when excluding the (mostly undetected) galaxies living 
inside Virgo. 

The value of M* for ETGs is ~ 5 times lower than the value 
typically found in large, blind H I surv eys (e.g.,[Zw aan et al. 2003 , 
Springob, Haynes & Giovanelli|2005| |Zwaan et al.|2005a||Martin| 
et al.|2010| ). Both our and these previous surveys are complete at 
such high M(H I) levels, so that comparing the M* values is correct. 
This confirms that ETGs contain smaller amounts of H I than other 
galaxies. 

A comparison between the low-mass end slope a of the ETG 
H I mass function and that found by previous authors for galax- 
ies of different morphology is complicated because of the different 
way samples are selected. Our sample is magnitude-limited while 
previous studies analyse samples of H l-detected galaxies with no 
selection on galaxy luminosity. This means that previous studies 
include fainter galaxies than we do. Such galaxies contribute to the 
low-mass end of the H I mass function, causing a to decrease. In- 
deed, they find a values of -1.2 to -1.4, i.e., the mass function 
keeps rising with decreasing M(H I) below M*, unlike in our ETG 



sample (see dotted line in Fig.jSj. We present a comparison between 
ETGs and spiral galaxies with the same luminosity selection in the 
next section. 



7 COMPARISON TO SPIRAL GALAXIES 

In Fig. [9] we compare the M(H I) and M(H i)/Lk distribution of 
ETGs to those of spiral galaxies. Red histograms refer to ETGs in 
the ATLAS^^ H I sample - solid line, dashed line and shaded his- 
togram have the same meaning as in Fig.jS] The blue solid line rep- 
resents the distribution of spirals belonging to the ATLAS^° parent 
sample, with the additional selection criterion 6 ^ 10 deg as for the 
ETGs studied here (Sec. [5]). Spirals in the ATLAS^^ parent sample 
have the same Mk selection as ATLAS^^ ETGs. 

For spiral galaxies the value of M(H I) is derived from the H I 
flux available in HyperLeda ( [Paturel et al.|2003| l, assuming the dis- 
tance given in Paper I. The Hi flux is available for 390/418 = 
93 percent of all spirals. From HyperLeda, it is not possible to 
know whether the remaining 28 spirals are undetected or unob- 
served in H I. However, about 80 percent of them are fainter than 
Mk= -22.5, so that they are likely to contribute to the low-M(H I) 
end of the spiral distribution. Furthermore, ~ 2/3 of them have 
distance larger than ~ 30 Mpc, which makes them difficult to de- 
tect depending on the depth of the data and may explain why they 
are missing from the HyperLeda database. This small fraction of 
missing galaxies and the relatively high uncertainty on HyperLeda 
H I fluxes are not particularly important for the type of comparison 
presented in this section. 

As discussed in Sec.[6] the M(H I) distribution of ETGs is very 
broad between (at least) a few times 10^ and a few times 10^ M©, 
while at lower H I masses it is uncertain because of many upper lim- 
its. Spirals in the ATLAS^^ parent sample look very different, with 
a much narrower distribution peaking at M(H l)~ 2 x 10^ M© and 
a tail towards H I masses below a few times 10^ M©. The number 
of galaxies in this tail might be slightly higher if most spiral galax- 
ies without available H I flux have small M(H I) (as argued above). 
However, this will not be a very large effect, and Fig. [9] shows that 
the relative number of galaxies with M(H I) below a few times 10^ 
M© is larger in ETGs than in spirals when the same Mk selection 
is applied. Therefore, for Mk< -21.5 the slope a of the spiral H I 
mass function might be much smaller than the value of ~ -0.7 de- 
rived for ETGs in the previous section. Overall, 86 percent of all 
ETGs and 87 percent of all spirals have M(H I) below and above 
5 X 10^ M©, respectively. 

The M(H i)/Lk ratio follows a peaked distribution in spirals 
owing to the known correlation between M(H I) and galaxy lumi- 
nosity ( |Shostak|I978| ). Such correlation does not hold for ETGs, so 
that their M(H i)/Lk distribution is much broader. In fact, it is ap- 
proximately flat between M(Hi)/Lk= 10~^ and 10~^ M©/L© and 
drops to zero above this value, where the distribution of spirals 
peaks. We find that 85 percent of all ETGs and 92 percent of all 
spirals have M(H i)/Lk below and above 10~^ M©/L©, respectively. 

Early-type galaxies and spirals are distributed very differently 
in Fig. |9] demonstrating that ETGs, as a family, contain much 
less neutral hydrogen gas than spirals (a similar indication comes 
from the low M* value of the best-fitting Schechter function de- 
scribed in the previous section). However, their distributions over- 
lap significantly in the range M(Hl)=5 X 10^ - 5 X 10^ M© and 
M(Hi)/Lk=3 X 10"^ - 10"^ M©/L©. A significant fraction of all 
ETGs contain as much H I spiral galaxies. 

This is in agreement with results of |Catinella et al.|(2010| ) as 
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Figure 9. Normalised histogram of M(H i) (left) and M(H i)/Lk (right) for ETGs and spirals in the ATLAS^^ parent sample. The dahsed red line refers to 
H i-detected ETGs. The solid red line is the corrected histogram (see Sec.[6j. The shaded histogram represents the distribution of upper limits for undetected 
galaxies. The solid blue line refers to spirals with available H I flux in HyperLeda. Error bars are obtained assuming a binomial distribution. 



part of the GASS survey of massive galaxies. They show a weak 
relation between M(Hl)/Mstar and the concentration index mea- 
sured from galaxies' optical images (their Fig. 8). This index is 
proxy for the bulge-to-disc ratio of a galaxy (e.g., | Weinmann et aT] 
|2009 X so that it is often used to select samples of early- or late-type 
systems. Adopting standard criteria for this selection, Catinella 
|et al.| ( |2010| ) results show that spiral galaxies are on average H l- 
richer than ETGs, and that most H I non-detections are ETGs for 
M(H i)iim/^star of a fcw percent. However, consistent with our re- 
sult, they also show that there is substantial overlap between the gas 
content of early- and late-type galaxies. 

Another element of similarity between spirals and ETGs with 
H I is that the majority of all H l-rich ETGs show disc-like H I mor- 
phology (64 percent over the entire sample and 80 percent above 
M(H l)=5 X 10^ Mo). Therefore, a significant population of ETGs 
exists which has similar M(H I), M(H i)/Lk and H I morphology as 
spiral galaxies. What is then the difference between the H I proper- 
ties of these two types of galaxies? 

To investigate this we compare the distribution of H I column 
density N(ll I) in detected ETGs and spirals. For the latter we use 
total-H I images constructed from WSRT datacubes as part of the 
WHISP survey ( [van der Hulst|2002] ). We make this comparison as 
fair as possible by: (i) using WHISP images at a 30 arcsec angular 
resolution, similar to that of our ETG H I images, and ( ii) studying 
only WHISP galaxies which either belong to the ATLAS^^ parent 
sample (54 galaxies), or fall within the same recessional velocity 
and Mb range as spirals in the ATLAS^^ parent sample but are 
outside the sky area covered by it (29 galaxies). 

We build the normalised N(ll I) histogram for each detected 
ETG and spiral separately. To do so we use for each galaxy all 
pixels with detected H I emission in the total H I image. Figure [T0| 
shows the median of all ETG histograms (red line) and of all spiral 
histograms (blue line). Both distributions increase with decreasing 
N(il I) down to the sensitivity limit of the observations. WHISP 
data cover a wide range in sensitivity and are on average a factor of 
~ 5 shallower than our observations of ETGs. The relative normali- 
sation of ETG and spiral N{ll I) distributions is therefore uncertain. 

We choose to normalise the histograms in the figure to the 
integral of the respective best-fitting Schechter function above 
N(tl l)= 10^^ cm~^. These best fits are shown by solid black lines 
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Figure 10. Distribution of H I column density N(ll i). The solid red line 
shows the median distribution for detected ETGs. The solid blue line shows 
the median distribution for spirals in the WHISP sample (see text). The 
shaded area corresponds to the 16th and 84th percentile (equivalent to 
±lcr for a Gaussian distribution). Thin black Hnes represent the best-fitting 
Schechter function for ETGs and spirals, with a dotted line showing the 
extrapolation of the best fit to column density below the sensitivity of the 
observations. The dashed red line represents the column density sensitivity 
of the shallowest total H i images presented in Fig. |Al| 



in the A/^(Hl) range where the fit was performed, and their ex- 
trapolation to lower column density is represented by a dotted 
black line. For spiral galaxies we find a = -0.99 + 0.04 and 
A^* = (1.03 + 0.07) X 10^^ cm-2. For ETGs A^* is an order of magni- 
tude lower, (9.2 + 0.6) x 10^^ cm~^. This value is close to our sen- 
sitivity threshold (red dashed line in the figure) so that a is not well 
constrained for ETGs. Therefore, we assume the same flat slope 
found for spiral galaxies (a ~ -1 J] We note that, although in dif- 



^ The best-fitting Schechter function to the A^(H i) distribution of spirals 
has a flat slope at low column densities not fully probed by the WHISP 
data. We verify the prediction of this fit at low column density by analysing 
a very deep WSRT H I image of a prototypical nearly face-on spiral galaxy, 
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ferent contexts, several authors have successfully parameterised the 
distribution of A^(H I) in galaxies using a Schechter function (e.g., 
IZwaan et al.|2005b| ). 

Figure[TO|and the above Schechter function parameters show a 
significant difference between ETGs and spirals. The N{Y{ I) distri- 
bution of spirals is very broad and stretches up to a few times 10^^ 
cm~^. In contrast, the distribution of A/^(Hl) in ETGs drops very 
quickly and, on average, stops at ~ 3x10^^ cm~^. Most ETGs never 
reach column densities above 5x10^^ cm~^ at the resolution of our 
data. This value is the mean column density of H I within the bright 
disc of spiral galaxies, derived from the tight relation between H I 
mass and Hi radius within the 1 Mo/pc^ = 1.25 x 10^^ cm~^ 
isophote (e.g., [Broeils & RheelpWTl [Noordermeer et al.|[2005] ). 
Therefore, ETGs rarely host neutral atomic gas as dense as even 
the average column density in the disc of spirals. It is interesting to 
note that the few ETGs with H I at about lO^^ cm-^ - NGC 2685, 
NGC 2764, NGC 3619 and NGC 7465 - are all galaxies with large 
amounts of molecular gas (Paper IV), complex and prominent dust 
distributions and star-forming regions (Paper II). 

So we find that, as far as Hi properties are concerned, the 
main difference between ETGs with large amounts of H I and spi- 
rals is that the former miss the high-column-density H I typical of 
the bright stellar disc of the latter. This is reasonable as the star- 
formation rate per unit area is much larger in spirals than in ETGs. 
Instead, the H I found in ETGs has column densities similar to those 
observed in the outer regions of spirals. In spirals, gas in these re- 
gions often exhibits warps and unsettled morphology and kinemat- 
ics, and is taken as a signature of the on-going accretion of gas on 
the host galaxy ( [Sancisi et al.|2008l ). This is similar to what we find 
in ETGs. As mentioned in Sec. [4] most of the H l-rich ETGs have 
a disc- or ring-like morphology, but in many cases part of the gas 
has not settled yet (or has been disturbed recently). This suggests 
that despite appearing as extremely different objects in their central 
regions, H l-rich ETGs and spirals may look very similar in their 
outskirts. 

Based on numerical simulations, several authors have argued 
that it is possible to form an ETG with a disc by merging two gas- 
rich galaxies (e.g., 'Na ab '&^Burk ert 2001, Barne s 2002, Springel 
& Hemquistl2 005 , Naa b, Jesseit & Burkert 2006| [Robertson et al. 
20061 )7 The distribution of the gas around the simulated merger rem- 
nants resembles in many aspects the H I systems discussed here. 
The size and mass of the distribution varies greatly depending on 
the properties of the merging galaxies and the merger geometry, and 
evolves with time since the merging. Both gas discs and rings can 
form as a result of these events, and gaseous tidal tails can extend to 
large distance from the remnant and be re-accreted at a later time. 
Our results suggest that, in order for the remnant to be an ETG, gas 
must settle on these discs in such a way to keep its column density 
low. Future comparison to simulations could explore this aspect in 
more details. 



whether these properties depend on galaxy luminosity and envi- 
ronment. 

Figure [TT] shows M(Hl) (top) and M(Hi)/Lk (bottom) of all 
galaxies in the ATLAS^^ H I sample as a function of galaxy mag- 
nitude Mk (left), large-scale environment density Eio (middle) and 
local environment density E3 (right). Values of Mk are listed in Pa- 
per I. Values of 2io and ^3 are listed in Paper VII (see Sec.|2]for the 
definition of these parameters). 

In the figure, triangles indicate upper limits on M(H I). Other 
markers represent different H I morphologies as explained in the 
legend on the top-right panel (see Sec. |4] for a definition of H I 
morphological classes). Red and blue markers represent galaxies 
inside and outside the Virgo cluster, respectively. 

We note that the ATLAS^^ H I sample shows no relation be- 
tween Mk and Eio or E3. Therefore, trends with environment can 
be separated from trends with galaxy luminosity. 

8.1 H I and galaxy luminosity 

There is no strong trend of H I properties with galaxy luminosity. 
Figure [TT] hints to a lack of H I at bright Mk. This is in agreement 
with results by |Catinella et al.| ( |2010| ), although a comparison of 
our results is not straightforward because of the different sample 
selections. In our case however, the small number of bright galaxies 
makes this trend relatively weak from a statistical point of view. 
For example, we consider the two sub-samples of galaxies above 
and below M(Hl)=5 x 10^ M©, respectively. According to a two- 
sample KS test, the hypothesis that their Mk distributions are drawn 
from a same parent distribution is rejected just at the ~ 93 percent 
confidence level. 

Massive galaxies appear to exhibit a disturbed H I morphology 
more frequently than fainter systems. The fraction of w's is 7/33 = 
21+8 percent for the brightest 20 percent of the ATLAS^° Hi 
sample, and 7/133 = 5 + 2 percent for the remaining galaxies. 
This seems to indicate that recent (major or minor) interaction with 
other gas-rich galaxies is relatively common for massive ETGs. For 
example. Paper IX discuss evidence of the major-merger origin of 
2 of these massive u galaxies (NGC 0680 and NGC 5557) based 
on the H I data presented here and on deep optical imaging, and 
argue that the formation event must have occurred after z ~ 0.5. 
However, while H I observations can reveal these past interactions, 
in a large fraction of these systems the H I is distributed at large 
distance from the galaxy so that it is not clear whether in all cases 
gas is or will be accreted on the central region of the ETG (e.g., 
NGC 3193 and NGC 5557 in Fig.|4}. 

Figure 8 suggests also that the fraction of J's is higher in 
fainter galaxies: 4/33 = 12 + 6 percent for the faintest 20 percent 
of the ATLAS^^ H I sample, against 6/133 = 5 + 2 percent for the 
rest of the galaxies. However, within the error bars this difference 
is only marginal. 



8 RELATION BETWEEN H I, GALAXY LUMINOSITY 
AND ENVIRONMENT 

In previous sections we describe the distribution of ETG H I mor- 
phology, M(Hl) and M(Hi)/Lk. In this section we investigate 

NGC 6946 jBoomsma et al.|2008) . We find that the A^(H i) distribution in 
this galaxy is flat and in excellent agreement with our best-fitting spiral 
Schechter function between a few times 10^^ and ~ 10^^ cm~^, confirming 
the validity of our fit. 



8.2 H I mass and environment 

Several previous authors reported a strong difference in H I detec- 
tion rate between ETGs living inside and outside the Virgo cluster. 
Within the ATLAS^^ Hi sample we detect 4/39 = 10 + 5 per- 
cent of all ETGs in Virgo and 49/127 = 39 + 6 percent of all 
ETGs outside Virgo. This is in good agreement with results from|di| 
[Serego Alighieri et al.H2007[ 2/32 = 6 + 4 percent inside Virgo for 
Mb < -18, comparable to our Mk selection), |Grossi et al.|p009[ 
5/14 = 36 + 16 percent outside Virgo with the same Mb selection) 
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Figure 11. M(H i) (top) and M(H i)/Lk (bottom) plotted against galaxy Mk (left), Eio (middle), and E3 (right). Different symbols represent different H I 
morphologies as explained in the legend. Red and blue markers represent galaxies inside and outside the Virgo cluster, respectively. 



and 010 on the Sauron sample (1/13 = 8 + 8 percent inside and 
14/20 = 70+ 19 percent outside Virgo; the latter is slightly higher 
than our result but the typical distance of galaxies in the Sauron 
sample is lower and, unilke OlO, our detection rate does not in- 
clude a number of cases where it is unclear whether H I belongs 
to the observed ETG or not - see Fig. |A2[ ). Our result represents a 
significant improvement on the statistical significance of these de- 
tection rates. 

While we confirm the Virgo vs. non- Virgo dichotomy, we note 
that environmental effects are more subtle than this. Firstly, M(H I) 
and M(H i)/Lk vary smoothly with environment density rather than 
showing clear breaks and sharp transitions. Figure [TT] suggests the 
presence of an envelope of decreasing M(H I) and M(H i)/Lk with 
increasing environment density over the entire density range cov- 
ered by our sample. 

One cause of this envelope is that the H l-richest ETGs live 
in the poorest environment. To investigate the significance of this 
result we compare the S 10 distributions of galaxies with M(H i)/Lk 
above and below 3 x 10~^ Mq/L©. A two-sample KS test rejects 
the hypothesis that the two distributions are the same at the 99.98 
percent confidence level (99.74 percent if we exclude galaxies in 
Virgo from this analysis). Similarly, the E3 distributions of the same 
two samples are different at the 99.6 percent confidence level. 

We find a similar result in our study of the CO content of 
ETGs (Paper IV) - the most CO-rich ETGs live in the poorest en- 
vironments. Therefore, the observation of both H I and H2 suggests 
that ETGs living in low-density environment accrete more cold gas 
(possibly from the surrounding medium) and/or can retain it for 
a longer period of time. This can be caused by the lack of both 
large neighbouring galaxies and a hot medium in such poor envi- 
ronments. 

Another clue to the different gas-accretion histories of ETGs 



living in different environments might be the way star forma- 
tion activity changes as a function of environment density. We di- 
vide the full ATLAS^^ sample in four Eio bins with bin edges at 
logio 2^10/Mpc^ = -2, -1,0, 1,2. Moving from the poorest to the 
richest environment we find dust/blue features (indicative of star 
formation; see Sqc.^, in 21/65, 18/105, 6/47 and 6/43 of all 
galaxies. We conclude that star formation signatures are more fre- 
quent at the poorest environment (32 + 7 percent) than in all the 
other Eio bins (~ 15 + 5 percent in each of them separately). 

Variations in the physical properties of the medium (e.g., den- 
sity and temperature) and in the rate of galaxy interaction might 
be causing the gradual decrease of M(H I) and M(H i)/Lk with en- 
vironment density. Such transition in H I properties is visible also 
within the Virgo cluster, suggesting that galaxies living at the out- 
skirts of Virgo are different from those closer to its centre. 



8.3 Centre and outskirts of the Virgo duster 

Galaxies detected in Virgo are NGC 4262, NGC 4406, NGC 4694 
and NGC 4710. Their projected distance from M87 is 1.0, 0.5, 1.3 
and 1.6 Mpc, respectively. Surface-brightness fluctuation distances 
are available for NGC 4262 and NGC 4406, and place them at a 
distance of 1.8 and 0.4 Mpc from M87 along the line of sight, re- 
spectively (Paper I). NGC 4406 is the only galaxy detected in H I 
close to the cluster centre, so that the H I detection rate in this re- 
gion is just 1/26 = 4 + 4 percent. In contrast, the detection rate 
outside 1 Mpc from M87 is 3/13 = 23 + 13 percent. So we find that 
H I detection rate, M(H I) and M(H i)/Lk (see Fig. [TT} in the cen- 
tre of Virgo are significantly lower than outside the cluster, while 
Virgo's outskirts appear as a transition region. 

In Paper VII we find that the variation of galaxy morpholog- 
ical mix with environment density also continues within Virgo. In 
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particular, the morphology-density relation changes slope at the 
Virgo cluster centre, where all Virgo slow rotators are found. This 
is in agreement with the above H I result and strengthens the con- 
clusion that different processes for the formation of ETGs must be 
at work inside the cluster core. These processes have the effect of 
generating a population of ETGs which is poorer of H I and signif- 
icantly richer of slow rotators than galaxies in less dense environ- 
ments, including the cluster outskirts. 

The above results are in good agreement with conclusions 
from the VIVA H I survey of spirals in Virgo by [Chung et al.| (|2009|). 
They find very H l-rich spirals with a gas disc larger than the op- 
tical disc (typical of field spirals) only further than 1 Mpc from 
M87 in projection. Closer to the cluster centre they find only galax- 
ies whose H I disc is truncated to the size of the stellar disc (or 
even smaller for systems very close to M87), and galaxies with gas 
tails caused by a combination of ram-pressure stripping and tidal 
interaction with neighbouring galaxies. The 1 Mpc cluster-centric 
radius around which [Chung et al. (2009 ) see a transition in the H I 
properties of spirals matches our result on the variation of ETG H I 
properties within Virgo. 

Neutral hydrogen is known to be quickly removed from galax- 
ies as they move within the hot intra-cluster medium (e.g., |Gunn &| 
Gott 1972 ). This is indeed the common interpretation of the fact 
that Virgo ETGs, which are a virialised population and have there- 
fore completed many orbits within the cluster, are so H l-poor (e.g., 
[di Serego Alighieri et al.|2007l OlO). Because of these considera- 
tions, the relatively high M(H I) and H I detection rate of ETGs in 
the outskirts of Virgo suggests that at least a fraction of them is 
falling now for the first time in the cluster (a similar argument is 
commonly used to claim that Virgo spirals are making their first 
passage through the cluster). This implies that at least a fraction 
of all Virgo galaxies may have joined the cluster already with an 
early-type morphology, rather than having become ETGs within the 
cluster. Their current morphology may result from pre-processing 
inside galaxy groups, as suggested also by our result that the main 
driver of galaxies' morphology-density relation are processes oc- 
curring at the galaxy-group scale (Paper VII). 

The H I properties of ETGs in Virgo are in contrast with results 
of our CO observations (Paper IV). Unlike H l-rich objects, ETGs 
with molecular gas do not avoid the cluster centre. Their detection 
rate is 5/33 = 15 + 7 percent within 1 Mpc from M87 and 6/25 = 
24+10 percent farther than 1 Mpc from M87 (for comparison, it 
is 47/202 = 23 ± 3 percent outside Virgo). Given the errorbars, it 
is still possible that the CO detection rate is lower in the centre of 
Virgo than at its outskirts or outside the cluster. However, it is clear 
that the difference between ETGs living in different environments 
is much more obvious in H I than in H2 (see jKenney & Young|1989| 
for a similar result on spiral galaxies). 

The H2 mass of ETGs inside the cluster is also consistent 
with that of galaxies outside it (with the aforementioned exception 
of the very H2-rich galaxies in the poorest environment; see Sec. 
[8.2[ ). This implies that the lack of H I in Virgo is related to a high 
M(H2)/M(H I) ratio rather than an overall gas deficiency. A possi- 
ble explanation is that H I is more easily stripped than H2 when a 
galaxy moves through a hot medium. This is in general true over 
the full extent of a galaxy because H I is typically distributed out to 
larger radii than CO and is therefore less bound to the galaxy ( [Fu-j 
[magalli et al.[2009[ ). This situation is particularly frequent in ETGs, 
where the H I often extends to tens of kpc from the galaxy (see Fig. 
[Al[ ). However, many ETGs outside the cluster contain H I in their 
central regions, so the question is why we do not detect such central 
H I (but do detect central CO) in Virgo. 



Neutral hydrogen in the centre of ETGs might be more easily 
stripped than CO because of its lower surface density ( Yamagami[ 
[& Fujita|[2011| ). Another important indication is that ETGs with 
CO appear as a virialised population within Virgo (Paper IV), so 
that they may have completed many orbits within the cluster. A 
prolonged residence within the cluster might be the reason why 
these galaxies are so H I poor. They are not only stripped of their 
gas, but the lack of other H l-rich companions means that no further 
Hi accretion can occur (see OlO). While this seems a plausible 
picture, we note that at least in some cases CO is stripped together 
with the H I from the centre of galaxies ( [Vollmer et al.|2008] ), so 
that we do not have a definitive answer to this problem yet. 

8.4 H I morphology-density relation 

In Sec. |8.2| we show that galaxies in poor environment can reach 
higher H I masses than those found at higher density. These ob- 
jects also tend to have a more regular H I morphology than gas-rich 
ETGs in denser environments. We illustrate the relation between 
H I morphological mix and environment density in Figs.[12|and[13| 
These show total H I images of the ten detections with the low- 
est and highest ^3, respectively, sorted according to increasing ^3 
(left to right, top to bottom). Figure pj] includes 3/4 of all Virgo 
detections, all 3 Ursa Major detections (NGC 3998, NGC 4026 and 
NGC 4111; H I in these galaxies is discussed also by jVerheijen &[ 
[ Zwaan[[2b01| and 4 galaxies living in relatively rich groups. The 
fraction of settled H I systems (i.e., D and d) is high at the lowest 
environment densities: 8/10 = 80 + 28 percent, with the only u sys- 
tem possibly being a misclassified d (see Sec.|4|. In contrast, it is 
very low at the highest densities: 2/10 = 20+ 14 percent. At high 
2^3 most H l-rich systems look disturbed (including some classified 
as D, like NGC 3998). 

The trend of Hi morphology with environment density is 
stronger using E3 than Siq. For example, galaxies with a very dis- 
turbed H I morphology figure among the ten objects with the lowest 
Eio (NGC 1023, NGC 7465); and the fraction of settled H I systems 
at the lowest Eio is the same as that at the highest Sio (5/10 against 
4/10). This indicates that processes influencing the H I morphology 
occur on a galaxy-group scale, and that a disturbed H I morphology 
is more clearly related to the presence of close neighbours than to 
an overall, large-scale environment over-density. 

It is interesting that a similar conclusion was reached in Paper 
VII based on the morphology-density relation of galaxies within 
the ATLAS^^ volume. There we found that processes occurring on 
a galaxy-group scale must be the main driver of the relation, which 
is tighter and steeper using S3 rather than Siq. The variation of H I 
morphology as a function of local environment density confirms 
the importance of such processes for the evolution of ETGs. 



9 SUMMARY AND CONCLUSIONS 

We study the H I properties of all 166 ETGs in the volume-limited 
ATLAS^^ sample above 6 - \0 deg and further than 15 arcmin 
from Virgo A. The sample includes galaxies within 42 Mpc and 
brighter than Mk=-21.5. This is the largest sample of ETGs with 
deep, interferometric H I observations to date. 

We confirm earlier findings that the H I detection rate of ETGs 
depends strongly on environment density ( di Serego Alighieri et al.| 
2007, Grossi et al. 2009 OlO). We detect 4/39 = 10 + 5 percent 
of all ETGs inside the Virgo galaxy cluster, and 49/127 = 39 + 6 
percent of all ETGs outside it. This is consistent with results from 
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Figure 12. The ten H I detections with the lowest E3. Images are sorted according to increasing S3, left to right, top to bottom. We refer to the caption of Fig 
|Al| for a description of the content of each image. In this figure the top-right scale bar indicates 10 kpc at the galaxy distance 



NGC 4026 iu) 



NGC 5198 (u) 



NGC 4111 (u) 



NGC 3998 (L>) 




NGC 4521 (c) 




NGC 4406 iu) 



2.0 X 10^^-2 



Figure 13. The ten H I detections with the highest S3 . Images are sorted according to increasing S3, left to right, top to bottom. We refer to the caption of Fig. 
|Al| for a description of the content of each image. In this figure the top-right scale bar indicates 10 kpc at the galaxy distance 



previous authors, although on a stronger statistical basis. The lower 
detection rate of ETGs in H I surveys of larger volumes is easily 
explained by their much poorer sensitivity. For examp le the M(H I) 
detection limit in the study of Toribio et"aD( 201ll is above 10^ 
Mq in most of the surveyed volume. Figure [9 shows that, at such 
sensitivity, their survey is expected to detect only ~ 10 percent of 
all ETGs. 

We classify galaxies according to their H I morphology. We 
find that 1/5 of all ETGs (and 1/4 of all ETGs outside Virgo) host 
H I distributed on a regularly rotating disc or ring. These systems 
represent the majority of the H I detections, 64 percent, so that if 
an ETG is detected in H I it will most likely host a gas disc or ring. 
Another 8 percent of all ETGs (26 percent of all detections) host H I 
in an unsettled configuration, while the remaining detected galaxies 



are surrounded by H I clouds scattered around the stellar body and 
not obviously associated to it. 

Although most detections can be easily classified within this 
scheme, we do find a number of galaxies intermediate (or in tran- 
sition) between classes. Indeed, we claim that the H I morphology 
of ETGs varies in a continuous way, going from very regular discs 
and rings to unsettled gas distributions, and from these to systems 
of scattered clouds which may or may not be associated to the host 
galaxy. This continuity may be related to the time passed since the 
last major episode of gas accretion or stripping relative to the H I 
orbital time in each galaxy. 

We divide regular, rotating H I systems in two classes based 
on their size relative to the stellar body of the host galaxy: D (large 
discs/rings - 24 galaxies) and d (small discs - 10 galaxies). The 
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former contain between 10^ and 10^^ M© of gas distributed out to 
tens of kpc from the stellar body. In half of these systems the H I 
is morphologically or kinematically misaligned with respect to the 
stellar body. In contrast, J's contain typically less than 10^ M© of 
H I confined within the stellar body and morphologically and/or 
kinematically aligned to it in nearly all cases. 

We investigate the role of H I in fuelling star formation within 
the host ETG. Confirming results by OlO, we find that the de- 
tection rate of signatures of star formation (molecular gas, dust 
discs/filaments and blue features) is high (~ 70 percent) for ETGs 
with H I within the stellar body. Namely, such features are found 
in allJ's, ~ 60 percent of al D's with central H I and ~ 50 percent 
of all unsettled H I systems with central H I. On the contrary, they 
are found in just ~ 15 percent of all ETGs with no central H I (or 
no H I at all at the sensitivity of our observations), highlighting the 
role of H I in enriching ETGs with material for star formation. 

The ISM in the centre of ETGs is dominated by molecular 
gas, whose mass is a factor of a few up to ~ 100 larger than the H I 
mass. Galaxies hosting a small H I disc {d) reach particularly high 
M(H l)/M(H2) values (larger than 10 in half of the cases), suggest- 
ing that they are forming molecular gas at efficiencies comparable 
to those of spirals. 

We parameterise the Hi mass function of ETGs with a 
Schechter function with M* = (1.8 + 0.7) X 10^ M© and a = 
-0.68 + 0.16. The value of M* is ~ 5 times lower than typical 
values found for samples dominated by spiral galaxies, confirming 
that ETGs are, as a a family, poorer of H I. The value of the faint- 
end slope a means that the mass function for a MK-limited sample 
decreases with decreasing M(H I) below M*. 

We compare the H I properties of ETGs to those of spirals. 
The M(H I) and M(H i)/Lk distributions of ETGs are very broad, 
reflecting the large variety of H I content of these galaxies, and con- 
firming the lack of correlation between ETG H I mass and luminos- 
ity. The distributions of spirals are much narrower and shifted to- 
wards larger M(H I) and M(H i)/Lk. Yet, we find a substantial over- 
lap in the M(H I) and M(H i)/Lk distributions of ETGs and spirals 
(consistent with |Catinella et al.pOlOj ). A significant fraction of all 
ETGs can have as much H I as spiral galaxies, and in the majority 
of the cases this is distributed on a large, rotating disc or ring. 

We investigate the difference between such H l-rich ETGs and 
spiral galaxies. We find that in the former the H I column density 
is typically very low. Early-type galaxies rarely reach the gas den- 
sity typical of the bright stellar disc of spiral galaxies, consistent 
with their lower star formation rate per unit area. The few galaxies 
that do, show clear signs of star formation and prominent dust fea- 
tures. Gas column density appears therefeore as the decisive factor 
in determining whether a galaxy is of early or late morphological 
type. 

We find a minor dependence of ETG H I properties on galaxy 
luminosity. Very luminous galaxies seem to contain less H I, but the 
only clear result is that H I in these systems is typically found in an 
unsettled configuration. This may reflect a higher rate of interaction 
with gas-rich companions for massive ETGs compared to fainter 
objects, although it is not clear how much of this gas is eventually 
accreted on the centre of the host galaxy. 

We do find clear trends of H I properties with environment. 
These go well beyond the known Virgo-vs.-non Virgo dichotomy 
mentioned at the beginning of this section. We find a smooth enve- 
lope of decreasing M(H I) and M(H i)/Lk with environment den- 
sity. Consistent with results from the ATLAS^^ CO survey pre- 
sented in Paper IV, we find that the gas-richest galaxies live in 
poor environments. This indicates that the cold-gas accretion rate is 



higher at these densities, and/or that these galaxies can retain their 
gas reservoir for longer periods of time. This effect is likely driven 
by a combination of a relatively quiet merging history and by the 
lack of a hot medium, and is accompanied by a higher fraction of 
objects exhibiting signs of recent star-formation in this environ- 
ment. 

The gradual decrease of H I mass with environment density 
continues within the Virgo cluster. Consistent with results from the 
study of spiral galaxies ( Chung et al.||200 9^, we find that the Hi 
properties of ETGs living at the outskirts of Virgo are intermediate 
between those of galaxies outside Virgo and those of the very H l- 
poor objects in its central region. Being H l-rich, galaxies at the 
Virgo outskirts are likely falling in the cluster for the first time, and 
they do so already with an early-type morphology. This suggests 
that pre-processing on a galaxy-group scale is fundamental for the 
evolution of ETGs, as concluded also in Paper VII. 

Galaxies inside ~ 1 Mpc from the centre of the cluster are 
significantly poorer of H I than galaxies outside Virgo, with a drop 
in the H I detection rate to just a few percent. This agrees with the 
observation of truncated H I discs in spiral galaxies in this region, 
and resonates with our finding that the fraction of slow rotating 
ETGs, which is very low at all environment densities, increases 
in the core of Virgo. This highlights that different processes must 
drive the evolution of galaxies deep inside the cluster, among which 
the presence of a hot medium (and therefore lack of H I). 

Unlike H I, molecular gas can be found in ETGs living deep 
inside Virgo. These systems have M(H2) similar to that of galax- 
ies outside Virgo, so that the lack of H I seems to imply a high 
molecular-to-atomic mass ratio. A possible interpretation is that 
molecular gas is more difficult to remove via ram-pressure strip- 
ping than H I, so that ETGs that have long been in the cluster (and 
therefore long been stripped of their H I) can still host molecular 
gas deep in the centre of their gravitational potential (Paper IV). 
However, it is not clear why H I is not detected in the centre of 
these galaxies, as observations show that CO can be stripped to- 
gether with H I in at least some cases ( |VollmeretaL 2008). 

We find that the Hi morphology also depends on environ- 
ment density. In poor environment the typical Hi detection ex- 
hibits a large, settled H I disc or ring. These systems must have 
been in place for many Gyr and indicate that very recent accre- 
tion/merging events have been minor. On the contrary, disturbed 
and unsettled H I morphologies are very frequent at higher envi- 
ronment density. Such H I morphology-density relation is clearer 
when measuring environment density on a local scale, indicating 
that processes occurring on a galaxy-group scale are driving factors 
for the evolution of ETGs. This agrees with the study of the kine- 
matical morphology-density relation presented in Paper VII, where 
we conclude that processes occurring on this scale determine the 
increase of the ETG fraction (and decrease of the spiral fraction) 
with environment density. 

As a concluding remark, we have demonstrated that a large 
fraction of ETGs contain significant amounts of H I. It is quite nor- 
mal for these galaxies to host neutral hydrogen gas, and H l-rich 
ETGs should not be regarded as odd or rare systems. On the con- 
trary, H I seems to be playing an important role in replenishing 
ETGs with cold gas and fuelling the formation of new stars and 
a stellar disc in these systems. The amount and morphology of this 
H I depend weakly on galaxy luminosity, one of the main parame- 
ters determining some of the most important galaxy properties such 
as their structure and stellar populations. Even trends with envi- 
ronment density, another fundamental parameter for galaxy evolu- 
tion, are characterised by a very large scatter. This scatter suggests 
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that the gas accretion history of ETGs varies widely from galaxy 
to galaxy even at fixed galaxy mass and environment density. A 
logical consequence of this large scatter is that ETGs are not such 
because they lack H I - many of them contain plenty of it. Instead 
we have shown that most ETGs must accrete gas in such a way that 
it remains at low column density - too low to support the higher 
levels of star formation seen in spiral galaxies. 
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Figure Al. Total-H i contours on top of DSS images of all H i-detected galaxies. Images are centred on the target galaxy and show a fixed area of 180 x 180 
kpc^ at the galaxy distance. Contour levels are N(R i)=No x 2" with Nq given at the bottom of each image and n = 0, 1, 2, ... Contours are coloured black to 
red, faint to bright. The beam is shown on the bottom left. In the cases where the ROl image reveals H I missing from the standard image we show its lowest 
contour with a thick grey line (contour level in parenthesis). 
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Figure Al. Continued 
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Figure Al. Continued 
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Figure Al. Continued 
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Figure Al. Continued 



© 2010 RAS, MNRAS 000,[T]|30l 



26 Paolo Serra et al. 

NGC 0770 




NGC 3226 















. ■ ^ I2.6 X 10i9cm-2 



NGC 3605 




NGC 4283 









■ 




• ■ . 








■■'.0 . . 


.■ 


■ 1.1 (0.9) X IQi^cm-^ 



NGC 2481 


NGC 2778 




+ ■ ■* ' 

' 0. . 

. " .■ . ^ ■ 

* ■ '■ ' * 


[0] . . 2.3 X IQi'^cm-^ 




_ . .■ 3.1 X lO^^cm-^- 


NGC 3379 


NGC 3599 


* 


. ■ ■ 

• 


. ■ ■ 2.0 X lO^^cm-^ 


[0] ■ • • 3.2 X 10i9cm-2- 


NGC 3607 


NGC 4168 


■ . * ■ 




[0] . 1.2 X IQi^cm-^ 




[0] ' - " 2.0 X 10i9cm-2 


NGC 5481 


NGC 6278 


. 


my . ^ 

*■ 


3.1 X 10^^-2 




.' , .2.4 X lO^^cm-^ 



Figure A2. Fields where H I is detected in the proximity of an ETG but is either clearly associated to a nearby galaxy or cannot be unambiguously associated 
to our target galaxy. See the caption of Fig. [AT] for a description of the images. 
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Table B 1 . Hi properties of ATLAS^^ galaxies. 



Name 


noise 


beam 


logio M(Hl) 


logio M(Hi)/Lk 


H I class 


notes 


source 




(mJy beain"^) 


(arcsec^) 


(Mo) 


(Mo/Lo) 








(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


(8) 


IC 0598 


0.46 


37.7X 37.3 


<7.45 


<c 


9 on 


- 


- 


- 


IC3631 


0.60 


77.8x31.8 


<7.71 




9 ZLO 


- 


- 


- 


NGC 0661 


0.51 


48.2X 32.9 


<7.37 


<c 


99 


- 


- 


- 


NGC 0680 


0.53 


51.4X 33.7 


9.47 




1 ^1 
— i .Z) i 


u 


- 


- 


NGC 0770 


0.55 


66.2X 33.5 


<7.56 


<c 


9 78 
— z. / o 


- 


- 


- 


NGC 0821 


0.30 


91.0x31.8 


< 6.91 


<c 


zL no 


- 


- 


s(4) 


NGC 1023 


0.20 


39.0X 34.6 


9.29 




1 


u 


- 


s(4) 


NGC 2481 


0.52 


49.1x 34.2 


<7.42 


<c 


— J.ZJ 


- 


- 


- 


NGC 2549 


0.43 


39.1x 36.4 


< 6.51 


<c 


— J. 1 o 


- 


- 


s 


NGC 2577 


0.48 


51.0X 34.0 


<7.35 


<c 


—J.J J 


- 


- 


- 


NGC 2592 


0.50 


46.8X 34.9 


< 7.18 


<c 


98 


- 


- 


- 


NGC 2594 


0.50 


46.8X 34.9 


8.91 




— 1 .J J 


D 


P,R 


- 


NGC 2679 


0.47 


41.5x36.3 


<7.35 


<c 


^ OS 


- 


- 


- 


NGC 2685 


0.23 


28.4X 25.6 


9.33 




— 1 10 


D 


W 


j(4) 


NGC 2764 


0.50 


52.3x 33.1 


9.28 




1 ^1 


D 


W,L 


- 


NGC 2768 


0.41 


37.5X 36.8 


7.81 




— J.JO 


u 


ROl 


s 


NGC 2778 


0.47 


40.2X 36.1 


<7.06 


<c 


"X IZL 
— J. IH- 


- 


- 


- 


NGC 2824 


0.49 


45.7X 34.4 


7.59 




9 8Q 


d 


ROl, A 


- 


NGC 2852 


0.47 


38.4X 36.3 


<7.27 


<c 


9 Q1 


- 


- 


- 


NGC 2859 


0.48 


40.5X 36.2 


8.46 




9 SO 


D 


R 


- 


NGC 2880 


0.48 


38.4X 38.3 


<7.03 


<c 


Al 

— J.H- / 


- 


- 


- 


NGC 2950 


0.47 


38.5x 37.1 


<6.69 


<c 


7Q 
—J. ly 


- 


- 


- 


NGC 3032 


0.24 


46.3X 35.0 


8.04 




9 08 


d 


C 


s(4) 


NGC 3073 


0.47 


38.6X 37.0 


8.56 




1 Af\ 

— 1 .H-D 


u 


- 


- 


NGC 3098 


0.50 


46.7X 33.5 


< 7.12 


<c 


^ 98 


- 


- 


- 


NGC 3182 


0.47 


38.3X 37.4 


6.92 




— J.U / 


d 


ROl 


- 


NGC 3193 


0.52 


51.0X 32.8 


8.19 




9 Q8 


u 


- 


- 


NGC 3226 


0.49 


54.8X 32.5 


<7.10 


<c 


—J.J i 


- 


- 


- 


NGC 3230 


0.62 


78.5X 32.6 


<7.71 


<c 


98 


- 


- 


- 


NGC 3245 


0.49 


42.0X 35.8 


<7.00 


<c 


% 7Q 
— J. 1 y 


- 


- 


- 


NGC 3248 


0.55 


49.6X 33.1 


<7.22 


<c 


-j.yj 1 


- 


- 


- 


NGC 3301 


0.53 


50.7X 32.8 


<7.13 


<c 


"X AQ 


- 


- 


- 


NGC 3377 


0.57 


74.8X 32.1 


<6.52 




— ^ 8Q 
— J. ay 


- 


- 


s 


NGC 3379 


0.59 


80.9X 32.1 


<6.49 






- 


- 


s 


NGC 3384 


0.59 


80.9X 32.1 


7.25 




"X Al 

— J.'-r 1 


c 


- 


s 


NGC 3400 


0.52 


42.2X 35.9 


<7.19 




—9 8S 


- 


- 


- 


NGC 3412 


0.60 


77.3X 32.2 


<6.55 


•< 


^ 78 
—J. 1 o 


- 


- 


- 


NGC 3414 


0.19 


44.9X 33.1 


8.28 




9 f\'X 


D 


/?01,P,R,L 


s(4) 


NGC 3457 


0.55 


61.8X 32.2 


8.07 




—9 00 


u 


- 


- 


NGC 3458 


0.48 


38.3X 37.2 


<7.35 




-3.21 


- 


- 


- 


NGC 3489 


0.25 


78.5X 35.0 


6.87 




-3.63 


d 


- 


s(4) 


NGC 3499 


0.50 


38.1x 36.9 


6.81 




-3.25 


d 


ROl 


- 


NGC 3522 


0.56 


58.0X 33.2 


8.47 




-1.51 


D 


/?01,P,R 




NGC 3530 


0.49 


38.4X 37.3 


<7.37 


< 


-2.74 








NGC 3595 


0.45 


39.5X 35.5 


<7.43 


< 


-3.19 








NGC 3599 


0.55 


61.0X 32.8 


<7.03 


< 


-3.17 








NGC 3605 


0.34 


62.9x31.0 


< 6.83 


< 


-3.21 






s(4) 


NGC 3607 


0.34 


62.9x31.0 


<6.92 


< 


-4.29 






s(4) 


NGC 3608 


0.34 


62.9x31.0 


7.16 




-3.61 


c 




s(4) 



© 2010 RAS, MNRAS 000,[T]|30l 



28 Paolo Serra et ah 



Table Bl (cont'd) 
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< /.Iz 


< — J.'^Z 






q; 


XT/^ /I /I Q ^ 

JNVjC 44 J j 


Q ^O 


1 on nw 1 on n 
loU.UXloU.U 


< /.ZJ 


< — J.DZ 






a 


NGC 4452 


4.65 


180.0x180.0 


< 7.27 


< -2.80 






a 


NGC 4458 


0.62 


39.4X 34.0 


< 6.91 


< -3.10 






s 


NGC 4459 


0.64 


69.7X 34.6 


< 6.91 


< -3.95 






s 


NGC 4461 


4.76 


180.0x180.0 


<7.33 


< -3.22 






a 


NGC 4473 


0.63 


36.7X 34.0 


< 6.86 


< -3.96 






s 


NGC 4474 


3.01 


180.0x180.0 


<7.08 


< -3.14 






a 


NGC 4477 


0.66 


71.5X 34.4 


< 6.95 


< -3.86 






s 
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Table Bl (cont'd) 



Name 


noise 


beam 


logio M(Hl) 


logio M(Hi)/Lk 


H I class 


notes 


source 




(mJy bearn"^) 


(arcsec^) 


(Mo) 


(Mo/Lo) 








(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


(8) 


M/^r^ A A CO 


n An 
U.4 / 


t^A 

D4.JX Jz.z 


< 0. /4 


< 


-3.21 








M/^r^ A Ad A 


U.J i 


/I C Q\/ C 


^ (1 Q/l 

< D.o4 


< 


-4.12 








Ksr^r^ A '^CM. 
iNLrL. 4jUJ 


'ill 


1 on Hn/ 1 on n 
ioU.UXioU.U 


^ n A A 
< /.i4 


< 


-3.46 






Of 


"VT/^ /I < o 1 


n /I Q 
U.4o 


QO QT A 

Jo.jX J /.4 


n n<. 




-3.13 


c 






M/^r^ A ^OQ 
INVjrC 4jZo 


T/l 

J. /4 


1 on 1 on n 
IoU.UXioU.U 


^ n A Q 
< /.io 


< 


-2.95 






Of 


INLrL. 4jjU 


U.DD 


70 n\/ QQ /I 


^ OQ 


< 


-3.33 






s 




C 1 A 

J. /4 


1 on n\/ 1 on n 
IoU.UXioU.U 


< /.jy 


< 


-2.80 






a 




U.oU 


01 /I \/ QO 1 
oi.4X Jz.i 


< D.O / 


< 


-4.16 






s 


iNLrL. 4jd4 


U.D / 


09 n\/ Q/i c 
oz. /X J4.J 


^ ^ Q1 


< 


-3.63 






s 


M/^r^ /I co^ 
JNLjL, 4j^^O 


J.Ui 


1 on nv/ 1 on n 
ioU.UXioU.U 


^ 'T 1 Q 
< /.iJ 


< 


-3.64 






a 


iNLrL. 4oUo 




1 on 1 on n 
ioU.UXioU.U 


< /.zz 


< 


-3.27 






a 


Nmr^ /I ^01 


U.D / 


1 7%/ Q/l Q 
oi. /X J4.J 


^ ^ 0^ 


< 


-4.11 






s 


INLrC 4djo 


Z.DJ 


1 on n>v/ 1 on n 
ioU.UXioU.U 


< /.iz 


< 


-3.39 






a 


M/^r^ Af^AQ 

INLrC 404^ 


Q HQ 


1 on n\/ 1 on n 
ioU.UXioU.U 


^ 'T 1 

< /.io 


< 


-4.32 






a 


iNLrL. 4DDU 


u.oy 


0/1 /I \/ Q/l Q 
o4.4X j4.j 


^ 00 
< D.OO 


< 


-3.51 






s 


INLtL. 4Dy4 


U.J / 


Q /I Tv/ Q 1 
j4. /X Ji.Z 


001 

o.zi 




-1.96 


u 




c 


IN LrC 4 / i U 


U.J / 


C\/ Q1 1 
d>'.jX JZ.i 


0/1 

D.o4 




-3.88 


a 






Mr^/^ /I '7'^'^ 
INLtL, 4 / JJ 


J. />' 


1 on ns/ 1 on n 
ioU.UXioU.U 


< /.iZ 


< 


-2.92 






Of 


JNVjL. 4 /j4 


J.JZ 


1 on n^/ 1 on n 
ioU.UXioU.U 


^ T 1 

< /.io 


< 


-3.59 






a 


iN LrC 4 / Dz 


j.Ui 


1 on (\\/ 1 on n 
ioU.UXioU.U 


< /.4U 


< 


-3.70 






a 


IN LrL. J i U J 


n /I Q 
U.4o 


Q7 /I \/ Q^ Q 
J /.4X JD.>' 


^7 
O.J / 




-1.69 


U 






JNLrC ji / J 


n /I ^ 
U.4d 


QO Q^ < 

Jo.jX jo. J 


Q QQ 
y.JJ 




-1.13 


U 






XT/^r^ ^ 1 QQ 


U.zj 


QT 1 \/ Q^ 1 
J /.iX Jj.i 


/IQ 




-2.46 


U 




s (2) 


INLrL. jz/j 


U.jU 


/in Q\/ Q/l Q 
4U.>'X j4.>' 


< D.oi 


< 


-3.45 










U.o / 


Q/l Q\/ Q/l n 

J4.JX j4.U 


< /.DJ 


< 


-3.34 






s 


iNLrC jJzZ 


n /I Q 
U.4o 


QO (\\/ QT ^ 
JO.UX J /.D 


< /.j4 


< 


-4.08 




A 

A 




iNLrL. JJ4Z 


U.J i 


QO n\/ QT /I 
JO.UX J /.4 


^ 7 ^n 
< /.jU 


< 


-2.85 










n An 
U.4 / 


QO n-^ Q^ < 
Jo. /X JO. J 


^ 7 /I < 
< /.4j 


< 


-3.90 




A 

A 






n /i'7 
U.4 / 


QO n\/ Q^ ^ 
JO. /X JO. J 


< / .JU 


< 


-2.77 










U.4 / 


QO n\/ Q^ ^ 
JO. /X JO. J 


< /.JZ 


< 


-2.60 










U.ji 


QO QO < 

jy.yx jy.j 


^ T Q/- 

< /.JO 


< 


-2.79 








iNLrC j4ZZ 


n An 
U.4 / 


QO Q\/ Q^ Q 
JO.jX JD.*^ 


7 07 




-2.92 


a 


7\Ui 




JNLtL. j4/J 


U.4o 


Q A^ Q^ 

Jo.4X JD.y 


< /.4U 


< 


-3.61 








JNLrL- j4/j 


n /IT 
U.4 / 


Q /I Q'7 n 

Jo.4X J /.U 


^ n 00 
< /.ZO 


< 


-3.19 








XT/^r"' ^/l Q 1 

INLiC j4oi 


n /IQ 


QO /I \/ QT n 

J0.4X J /.U 


< /.zi 


< 


-3.18 








XT/^ < /I o < 


n /IT 
U.4 / 


QT 0\/ QT < 

J /.oX J /.J 


^ n ^n 
< 1 A 1 


< 


-3.59 








XT/^r^ <<rvn 


(\ A(^ 

U.4d 


QT Ov/ Q^ 


^ n Q/r 
< /.JO 


< 


-2.73 








XT/^r^ ^^^T 

INLtL, jjj / 


n An 
U.4 / 


QQ 0\/ Q^ Q 

jy.zX jj.y 


^7 
O.J / 




-2.69 


u 






XT/^/^ <<0O 

INLrC JJoZ 


n /1 
U.4o 


QO Q^ 1 
Jo.jX JD.i 


Q ^< 




-0.97 


u 


T) 

K 






n /I Q 
U.4o 


/in n\/ Q^ ^ 
4U.UX jj.D 


< /.ij 


< 


-3.04 








XT/^r^ ^^Qi 


n /I Q 
U.4o 


Q'7 0\/ Q^ A 
J /.ZX JO. 4 


OQ 

o.oy 




-1.90 


T\ 

U 






XT/^/^ ^A.Qn 
JNLtL. jdo / 


U.JO 


QO A^ Q^ 


^ T QO 
< /.JZ 


< 


-3.28 








NCjC 5866 


0.47 


00 n 1 

38. jX 37.1 


6.96 




-3.95 


d 






NGC 6149 


0.53 


56.0X 32.4 


<7.56 


< 


-2.79 








NGC 6278 


0.52 


49.2X 32.9 


<7.67 


< 


-3.32 








NGC 6547 


0.52 


45.9X 33.3 


<7.63 


< 


-3.12 








NGC 6548 


0.53 


58.6X 32.3 


<7.12 


< 


-3.47 








NGC 6703 


0.46 


38.3X 36.4 


<7.18 


< 


-3.67 








NGC 6798 


0.60 


38.1x 37.0 


9.38 




-1.34 


D 


C 
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Table Bl (cont'd) 



Name 


noise 


beam 


logio M(Hl) 


logio M(Hi)/Lk 


H I class 


notes 


source 




(mJy bearn"^) 


(arcsec^) 


(Mo) 


(Mo/Lo) 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


XT/^/^ TO on 


U.OU 


^/l T 

d4.jX jz. / 


T OO 


— Z.JZ 


u 


Dm 

K\JL 




rsvjL. IddZ 


Ci /in 


/I 1 V/ Q 1 fk 

40. iX Ji.U 


D.DZ 


/I 1 Q 


c 




s (1.^) 


iNLrL. /4j4 


U.jj 


DJ.iX JZ.J 


^ T 1 ^ 


< — j.jj 








XT/^ T /I ^ T 

JNLrC /4j / 


f\ /lO 

U.4y 


J0.4X JJ.J 


< D.Oi 


< —J. DO 






s 


INLrC /4dj 


U. j4 


/T/r /Tv/ Q 


Q QO 


n /i^ 
— U.4d 


u 






rvji^ Uzooo / 


n ^T 
U.J / 


QC /I w ao 1 
OJ.4X jZ. i 


T 


9 CT 
—Z.J / 


c 








U.JO 


TO Q\y Q 1 T 

/o.oX Ji. / 


^ T /TO 

< /.DO 


< — Z.JU 




A 

A 




rvjC UJ J / j4 


U.ji 


/in q 
4U.jX jj.'s' 


< /.JO 


^ O /IQ 








rVjL. U444JJ 


U.Jo 


T^ /I w a 1 Q 


< /.DO 


^ 1 cc 
< — Z.JJ 










n /IT 
U.4/ 


QT Ow QT < 

J /.oX J /.J 


^ T CI 
< /.Ji 


< —Z.J / 








rOC Uj i / J J 


n /I ^ 
U.4j 


QQ Cw O/r 
JO.jX JD.D 


< /.JZ 


< — Z.JD 








rVjrl^ UDi4Do 


n ^/i 
U.j4 


^T 1 \/ '^0 A 
J /. iX jZ.4 


^ T C/l 

< /.j4 


^ A^ 
< — Z.4j 








PGC 071531 


0.52 


54.7X 32.5 


< 7.37 


< -2.64 








UGC 03960 


0.46 


50.0X 33.6 


7.79 


-2.28 


D 


W,L 




UGC 04551 


0.47 


39.7X 38.4 


< 7.25 


< -3.23 








UGC 05408 


0.47 


38.3x 37.1 


8.52 


-1.60 


d 






UGC 06176 


0.52 


51.6X 33.3 


9.02 


-1.36 


D 


W 




UGC 08876 


0.48 


39.4X 35.7 


< 7.43 


< -2.83 








UGC 09519 


0.48 


40.2X 35.6 


9.27 


-0.83 


D 


M 





Note. — 

Column (1): Principal designation from LEDA. 
Column (2): Median noise in the H I cube. 
Column (3): Major and minor axis of the beam. 

Column (4): Total H I mass calculated assuming galaxy distances given in Paper I. 
Column (5): Ratio between M(H I) and the ^-band luminosity given in Paper I. 
Column (6): H I class: D = large disc/ring; d = small disc/ring; u = unsettled; c = cloud. 

Column (7): Notes on H I data, morphology and kinematics: ROl indicates that additional Hi is detected in the 
ROl cube, and the M(H I) value is obtained from it; A = H I detected (also) in absorption; C = H I counter-rotating 
relative to the stellar kinematics; L = lopsided Hi morphology; M = Hi misaligned relative to the stellar 
kinematics; P = polar; R = ring; W = warp. 

Column (8): Source of H I data if other than the present study: c = [Chung et al.| ( |2009] ); j = |J6zsa et al.| ( |2009| ; o = 
Oosterloo; s = Sauron (M06, OlO); a = ALFALFA. The number in parenthesis indicates the integration time (in 
units of 12 h) for data obtained with WSRT observations longer than 12 h. 
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